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The use of CO. Laser technology for spacecraft tracking applica- 
tions is motivated by: 

(a) the high ranging accuracy potential associated with the 
large carrier and modulation frequencies involved, 

(b) its compatibility with CW analog modem methods that permit 
independent, range and range rate data extraction in a unified 


DD , 1473 EDITION Ol I NOV &S IS GbbOLETt 


UNCLASSIFIED 

SECURITY CLASSIFICATION OF IMIS PACE (W>»n Vaia £<ifara«/j 











1 I 4 II 111 


UNCLASSIFIED 

btCuHlTt LI AVjlMLAtlON Of Tm^> PALtlHI.wi Uulu hnlviuJi 

system concept. 

This report evaluates the conceptual design and error performance 
of a C0 2 laser ranging system on this basis. The study identi- 
fies ranging signal and subsystem processing alternatives, and 
their comprehensive evaluation yields preferred candidate solu- 
tions which are then analyzed to derive range and range rate 
error contributions. The performance results are presented in 
the form of extensive tables and figures that identify the 
ranging accuracy compromises as a function of the key system 
design parameters and subsystem performance indeces. The ranging 
(errors obtained are noted to be within the high accuracy require- 
ments of existing NASA/GSFC missions with a proper system design. 
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SUMMARY 

This report presents the results of a C0 2 Laser Ranging Systems 
Study. The use of laser ranging techniques is motivated by the 
high-accuracy satellite tracking specifications (0.05 meters in range, 
0.005 cm/sec in range rate) required in various NASA/GSFC missions. 
The use of CO 2 laser technology features high carrier (2.8 x 10 1 Hz) 
and modulation frequencies for accurate ranging, plus CW analog 
modem compatibility for a unified independent range and range rate 
data extraction. The range data can be obtained from the propaga- 
tion delays or phase offsets exhibited by a ranging event associated 
with the modulation signal, while the range rate data can be ob- 
tained from the doppler effect present in the received carrier 
signal. The independent extraction of range and range data has 
proven to be a desirable feature in tracking systems from orbital 
determination considerations. 

The use of harmonic sidetone or pseudo-noise code ranging 
principles have been effectively employed in independent range and 
range rate tracking systems with the ranging event being identified 
as the tone zero-crossing or the code epoch. The tone signal in- 
volves a simpler subsystem processing for acquisition and tracking, 
though a set of harmonic tones is required to assure both fine 
accuracy and ambiguity resolution in the range data. The code signal 
provides both these features in a single signal, but at the expense 
of acquisition time and code detectability compromises. The use of 
a hybrid signal employing a high frequency tone for fine ranging 
purposes along with a reduced set of lower tones plus one code for 
ambiguity resolution optimizes the ranging signal design rationale. 
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The high accuracy tone and the ambiguity resolving signals (lower 
tones, code subcarrier) should modulate the carrier in a sequential 
mode, so as to minimize the power sharing compromise between the 
carrier component (range rate signal) and the modulation sidebands 
(range signal) in the ranging signal. 

The choice of analog modem offers the conventional modulation 
options (amplitude and/or angle modulation, double or single side- 
band) , plus coupling modulation as an available candidate for laser 
ranging systems. The receiver processing also offers the options 
of incoherent envelope detection or coherent product demodulation 
for the baseband signal extraction. The use of narrowband phase 
modulation along with coherent product demodulation are shown to 
provide an effective carrier vs. sideband power sharing compromise, 
as well as an efficient: nonlinear processing that also exploits the 
available coherent carr.er at the receiver. The use of double side- 
band signaling is also noted to provide error propagation protection 
when compared to a single sideband signal alternative. The PM 
carrier vs. sideband power sharing is quantitatively illustrated as 
a function of the modulation parameters. It is noted that a maximum 
phase modulation index of 1 radian assigns a minimum 58.557, of the 
available power to the carrier component and a maximum 38.737, to 
the major tone sidebands, but an extension to 1.4 radian would 
reverse the power sharing to a 32,147, minimum for the carrier and 
a 58.737, maximum for the one sidebands. Such extended modulation 
index range is maintained in our C0 2 laser modulation evaluation, 
since it provides a signi icant improvement in the range vs. range 
rate accuracy compromise nduced by the sideband vs. carrier power 
sharing . 
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The carrier and tone signals are tracked at the receiver via 
phase-locked loops. These loops involve a design compromise based 
on the joint interest in wideband tracking for effective doppler 
dynamic and oscillator jitter reproduction vs. narrowband tracking 
for additive thermal noise rejection. This tracking compromise is 
evaluated for 2nd and 3rd order loop realizations involving practical 
non-ideal integrators in their loop filters, and their loop design 
rationale and locking error performance are established. The range 
and range rate error contributions of these locking errors are then 
formulated after modeling the data extraction operations performed. 
The error analysis thus permits the tracking loop design to be 
guided by either a threshold locking error rationale or a tracking 
data error rationale. In either case, the loop design goals are to 
provide accurate range and range rate data, while maintaining the 
locking error below established threshold bounds to assure the signal 
reproduction . 

The carrier loop design and range rate error performance is 
evaluated to illustrate that 3rd order loop realizations are required 
to maintain lock for the SNR conditions and orbital dynamic magni- 
tudes expected in the missions. The additive thermal noise error 
contribution to range rate data is shown to be rather small, with 
oscillator noise effects yet to be evaluated. The tone loop design 
and range error performance is also evaluated assuming two possible 
design approaches. The first case chooses a dynamic locking error 
bound ( y m degrees) to assure maximum doppler rate lock, and evaluates 
the corresponding dynamic ranging error bound (y, meters) as a func- 
tion of the major tone frequency (F m MHz). The second case reverses 
the design by choosing (p r ) and evaluating (w„) for each (F m ) . In 
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both cases, the rms locking error (o„ degrees) and rms ranging 
error ( o, meters) due to additive noise are evaluated as a func- 
tion of the phase modulation index (3 radians) and the major tone 
frequency (F m MHz) . The use of either 2nd or 3rd order loco reali- 
zations maintains lock for the tone loop design, and both design 
cases above are evaluated for both realizations in question. The 
results are tabulated and plotted a* families of curves that show 
(a r ) vs. (6) as a function of (F m ) , with each family corresponding 
to one design approach (y m or y, bound) and one loop realization 
(2nd or 3rd order). 

The tables that follow present some of the additive (thermal) 
noise range error magnitudes obtained under various logical design 
conditions, with oscillator noise effects yet to be evaluated. The 
tabulated results illustrate that small rms range errors within the 
accuracy specifications can be obtained by proper selection of the 
modulation frequency, modulation index, and loop order. The rms 
range error improvement capabilities of 3rd-order loop realizations 
are evident in each case. The functional dependence between the 
peak dynamic locking and ranging errors (y m , y r ) is indicated at the 
end of each table to permit verification of the loop locking status 
and the dynamic range error contribution in each case. 
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RMS Additive-Noise Range Errors in Centimeters: - 1 deg. 


(a) 

Case of y m ■ 1 deg 

and 2nd-order loop. 


Fm_ 

0=0.5 

0=1.0 

0=1.5 

1 

27.62 

15.21 

11.99 

3 

12.12 

6.670 

5.261 

10 

4.911 

2.704 

2.133 

30 

2.155 

1.186 

0.9355 

100 

0.8734 

0.4808 

0.3792 

300 

0.3831 

0.2109 

0.1664 


(b) 

Case 

of y m = 1 deg and 

3rd-order loop. 


Fm_ 


0=0.5 

0=1.0 

0=1.5 

1 


9.618 

5.295 

4.176 

3 


4.219 

2.323 

1.832 

10 


1.710 

0.9416 

0.7426 

30 


0.7503 

0.4131 

0.3258 

100 


0.3041 

0.1674 

0.1321 

300 


0.1334 

0.07345 

0.05793 

Note : 


100 

= » peak 

2.4 F m 

dynamic ranging error 

in cm. 
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RMS Additlve-Noiae Range Errors in Centimeters: yj, - 5 deg. 


(a) Case of y B ■ 5 deg and 2nd-order loop. 


fk. 

B=0.5 

8-1.0 

8=1.5 

1 

18.47 

10.17 

8.020 

3 

8.102 

4.461 

3.518 

10 

3.284 

1.808 

1.426 

30 

1.441 

0.7932 

0.6256 

100 

0.5840 

0.3215 

0.2536 

300 

0.2562 

0.1411 

0.1113 


(b) 

Case of yi m = 5 deg and 

3rd-order loop. 


Fm_ 

8=0.5 

8=1.0 

6=1.5 

1 

6.432 

3.541 

2.793 

3 

2.822 

1.553 

1.225 

10 

1.144 

0.6297 

0.4966 

30 

0.5017 

0.2762 

0.2179 

100 

0.2034 

0.1120 

0.08832 

300 

0.08922 

0.04912 

0.03874 


100 

Note: u r = = peak dynamic ranging error in cm. 

2.4 F m 


RMS Additive-Noise Range Errors in Centimeters: y r 


1 cm. 


(a) Case of y, * 1 cm and 2nd-order loop. 



6*0.5 

6-1.0 

6-1.5 

l 

70.16 

38.62 

30.46 

3 

23.39 

12.87 

10.15 

10 

7.016 

3.862 

3.046 

30 

2.339 

1.287 

1.015 

100 

0.7016 

0.3862 

0.3046 

300 

0.2339 

0.1287 

0.1015 


(b) 

Case of p r =1 cm 

and 3rd-order loop. 


5u 

6-0.5 

6-1.0 

6*1.5 

l 

24.43 

13.45 

10.61 

3 

8.144 

4.483 

3.536 

10 

2.443 

1.345 

1.061 

30 

0.8144 

0.4483 

0.3536 

100 

0.2443 

0.1345 

0.1061 

300 

0.08144 

0.04483 

0.03536 

Note 

2.4 y, F. 

Vm = 

100 

- peak dynamic locking 

error in deg. 
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RMS Additive-Noise Range Errors in Centimeters: y f 


5 cm. 


(a) Case of y r "5 cm and 2na-order loop. 


ZL 

6-0.5 

6-1.0 

6-1.5 

l 

46.92 

25.83 

20.37 

3 

15.64 

8.610 

6.791 

10 

4.692 

2.583 

2.037 

30 

1.564 

0 8610 

0.6791 

100 

0.4692 

0.2583 

0.2037 

300 

0.1564 

0.08610 

0.06791 


(b) 

Case of u, “5 cm and 

3rd-order loop. 



6-0.5 

6-1.0 

^1-5 

l 

16.34 

8.995 

7.094 

3 

5.446 

2.998 

2.365 

10 

1.634 

0.8995 

0.7094 

30 

0.5446 

0.2998 

0.2365 

100 

0.1634 

0.08995 

0.07094 

300 

0.05446 

0.02998 

0.02365 


2.4 y, F m 

Note: u m - ■ peak dynamic locking error in deg. 

100 
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General Introduction 

The use of CO 2 laser ranging techniques for high-accuracy 
satellite tracking applications is presently under consideration 
at NASA/GSFC. Various existing missions involving the GEOS and 
TIMATION spacecrafts are characterized by rather tight tracking 
accuracy requirements, such as 0.05 meters in range and 0.005 cm/ sec 
in range rate. These magnitudes represent an improvement in the 
state-of-the-art ranging capabilities when compared to existing 
tracking systems. The use of C0 2 laser technology for these pur- 
poses represents a natural solution approach, since it is compati- 
ble with various analog modem methods, so that one can both exploit 
the ranging accuracy advantages of high modulating frequencies while 
maintaining a CW carrier component at 2.8 x 10 13 Hz for accurate 
independent range rate data extraction. The availability of inde- 
pendent range and range rate data for orbital determination purposes 
has proven to be a desirable feature in tracking system design from 
error performance considerations. 

A study program of the C0 2 laser ranging system was conducted 
to this effect. The objective was to perform a comprehensive evalua- 
tion of the available modulation methods and tracking system design 
for C0 2 laser ranging, and then to establish their fundamental 
ranging accuracy limitations when specific orbital dynamics to be 
tracked and representative SNR level conditions are assumed. The 
study thus includes both a qualitative analysis to establish the 
relative merits of various analog modems and subsystem processing 
alternatives, as well as a quantitative evaluation of the tracking 
and ranging error performance of the system. 
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The analytical evaluation and error performance results of 
this study program are presented in this report. The ranging system 
design principles that guide and govern the conceptual design and 
development of high-accuracy tracking systems are first presented 
to establish a natural foundation. The characteristics of indepen- 
dent range and range rate data extraction are identified, and the 
ranging signal design and subsystem processing considerations that 
follow are discussed for our laser ranging application. The range 
vs. range rate accuracy compromise involved is emphasized to moti- 
vate a careful attention to the sideband vs. carrier power sharing 
compromise in the analog modem candidates. 

The ranging signal structure composition is analyzed to evaluate 
available options, and account for the fine ranging vs. ambiguity 
resolution operations involved. A hybrid signal composition that 
features the relative advantages of both harmonic sidetone and 
pseudo-noise code ranging principles is noted to optimize the ranging 
signal design compromises, along with the sequential application of 
the fine ranging and ambiguity resolution operations. The analog 
modem alternatives are carefully analyzed to establish their power 
sharing compromises and processing limitations that affect their 
ranging accuracy performance. The discussion includes all conven- 
tional AM, PM and AM/PM approaches, plus coupling modulation as an 
added option characteristic of C0 2 lasers. The power sharing and 
ranging error implications of envelope (incoherent) vs. product 
(coherent) demodulation methods are also evaluated for the modulation 
approaches in question. The specific limitations added by the laser 
modulation parameter values achievable are included in the discussion. 
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The tracking performance evaluation is specified for phase- 
locked extraction of the carrier (range rate) and tone (range) 
components of the ranging signal. The tone design parameters, loop 
design compromises, and tracking error analysis are identified for 
2nd and 3rd-order loop realizations. The analysis permits a clear 
comparison of 2nd vs. 3rd-order loop performance, plus emphasizes 
the loop design and error performance distinctions between theoreti- 
cal and practical limitations. The 3rd-order analysis features 
state-of-the-art technology, plus a unified formulation that readily 
guides natural extensions and potential improvements. The performance 
analysis is concentrated on orbital dynamics tracking in the presence 
of additive thermal noise, with the acknowledgement that oscillator 
noise effects should be ultimately included and may be significant 
with the high modulation and laser carrier frequencies involved. 

The range and range rate data errors contributed by the loop 
tracking errors are next evaluated. The analysis requires the care- 
ful modeling of the range and range rate extractor processing to 
properly reflect the effect of additive phase disturbances on the 
raw data samples obtained. In particular, two range rate extractor 
realizations are detailed and analyzed to derive the corresponding 
data error expressions. The report concludes with a quantitative 
error analysis and discussion based on these results, using typical 
link parameter specifications to derive the numerical error values. 

Two different loop design approaches are exhaustively evaluated to 
account for both locking and ranging error bound options in the 
formulation. All the numerical results are tabulated for reference 
purposes, plus families of characteristic performance curves 
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are included to facilitate the error analysis and logical design of 
the C0 2 laser ranging system. A glossary of the basic symbols used 
in this report is included at the end of the report for reference 
purposes and easy accessibility. 
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2 . Ranging System Design Principles 

The conceptual design and development of a spacecraft tracking 
system is based on the measurement of certain parameters that re- 
flect changes in the spacecraft trajectory, along with an orbital 
determination program that processes the measured data after some 
preliminary smoothing or editing operations. A conventional set 
cf measured parameters employed at NASA/GSFC for high-accuracy 
tracking purposes consists of the slant range and/or range rate of 
the spacecraft relative to a ground station network used as refer- 
ence points. The number of stations involved depends on the number 
of simultaneous measurements to be performed, the orbital dynamics 
to be tracked, and the accuracy performance to be obtained. 

The ground-centralized tracking systems extract the tracking 
data at the ground stations from the ranging signals received from 
the spacecraft. These signals may have been originally transmitted 
by the stations and returned by the spacecraft (two-way tracking 
systems) or they may originate at the spacecraft itself (one-way 
tracking systems) . Our application specifies a two-way ground- 
centralized tracking system, where the spacecraft returns a 2.8 X 
10 1 3 Hz laser ranging signal back to ground by means of a passive 
reflector. The simplified signal turnaround that results from a 
passive transponder bypasses the signal acquisition and processing 
issues that arise with an active transponder. 

The range data (r) may be obtained from the round-trip propaga- 
tion delay (t p ) of a suitable ranging event associated with the 
ranging signal. The time elapsed between the transmission and re- 
ception of the ranging event is measured at ground, and the range 
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c 

is derived via the relation r - 2 (t*). where (c) is the speed of 
light. The repetition period of the ranging event must be obviously 
smaller than the round-trip propagation delay to provide unambig- 
uous consecutive range measurements, in the absence of other avail- 
able information. However, the ranging signal may include ambiguity- 
resolving features to relax this requirement. In any case, epoch 
variations in the form of signal processing delays or added timing 
disturbances will contribute a range measurement error that must 
be calibrated out if possible or acknowledged as an accuracy limi- 
tation otherwise. Any departure from free-space propagation in the 
space channel will also affect the conversion of the measured delay 
into true slant range data. 

The range rate data (r) may, in turn, be obtained from the 
round-trip doppler effect (f d ) exhibited by a ranging signal fre- 

# p 

quency (f), via the relation r = f.. The use of a large frequency 

2f 

value is motivated for the attenuation of various error contributions 
in the doppler measurement process, and the link, frequency assign- 
ment in our application is indeed compatible with this spirit. The 
presence of frequency drifts will, in general, contribute an apparent 
doppler effect, and induce a range rate measurement error unless 
compensated for in the data extraction method. The doppler measure- 
ment itself may be affected by the presence of random phase jitter 
disturbances accompanying the ideal signal phase fluctuations, and 
unavoidable range rate errors will be introduced as a consequence. 

A departure from free-space propagation in the space channel will 
also serve as an accuracy limitation in the conversion of the measured 
doppler effect into true slant range rate data. 
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The derivation of range or range rate data from one another is 
obviously possible. However, the association of the aforesaid 
ranging event and doppler effect with different components in the 
ranging signal structure can provide range and range rate data 
independently, which often represents a desirable feature from 
error control and information reliability considerations. An un- 
avoidable compromise in the power allocation and accuracy perfor- 
mance then results. The available signal power is distributed 
among the range and range rate signal components, so that their 
respective receiver SNR's can only be improved at the expense of 
one another. The net effect is that the additive noise contribu- 
tion to range error decreases while its range rate error contribu- 
tion increases, and vice versa, as various power sharing percentages 
are considered. 

The problem of specifying quantitative preferences to optimize 
the independent range vs. range rate accuracy compromise becomes 
rather complex in general. One would have to account for the various 
options available in orbital dynamics, data rates, error magnitudes, 
orbital determination procedures, and state-of-the-art technology. 

The usual approach consists of considering each application in ques- 
tion, acknowledging the existing tracking accuracy capabilities as 
well as the projected improvements, and then investigating the 
variation of the range and range rate error magnitudes as a function 
of the subsystem design parameters to be specified. Such approach 
can indeed identify diminishing return conditions, where an incre- 
mental error improvement in one data is accompanied by a severe 
error degradation in the other. The method then yields a range of 
design parameter values that represent favorable choices, though 
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the order of magnitude of the range and range rate error requirements 
is assumed a-priori to guide the compromise evaluation. 

The conceptual design of the ranging system thus begins with 
the identification of the ranging signal event and the carrier fre- 
quency selected to respectively provide the independent range and 
range rate tracking data. The ranging signal design then considers 
the integration of these features into a composite signal propaga- 
ting through the space link. The ranging signal structure is 
specified following a careful evaluation of modulation methods 
available to the system designer, which concludes with a choice 
among the preferred candidate solutions and the assignment of nu- 
merical values to the modulation parameters involved. The parameter 
assignment involves an unavoidable compromise between the range 
and range rate accuracy performance of the system, since the avail- 
able power budget must be shared by the ranging signal components 
providing the independent range and range rate data. 

The use of sidetone zero-crossings or pseudo-noise code for- 
mats as the ranging events is compatible with the high- accuracy 
long-range requirements of our laser tracking application. The 
doppler effects present in the carrier or clock components associated 
with these signals are also capable of providing a high range rate 
resolution. These ranging principles have been reliably employed 
in the design and development of tracking systems for space applica- 
tions, and represent state-of-the-art procedures for CW ranging 
purposes. These procedures are indeed compatible with C0 2 laser 
modulation principles, which can then benefit from the technological 
experience and facilities already available. 
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The selection of a ranging signal structure is also linked to 
the problem of subsystem specification, since the ground statio.. 
must generate, transmit, receive and track the ranging signal, plus 
extract and process the range and range rate data. The ground 
oscillators used in the transmitter signal generation and receiver 
signal processing contribute long-term frequency drifts and short- 
term instabilities that ultimately limit the tracking performance 
and accuracy capabilities of the system. The analysis of oscillator 
noise effects is, of course, emphasized by the high frequencies 
(carrier and modulation) involved in our laser ranging system. The 
signal processing itself is equally critical in the wideband applica- 
tion envisioned, since the design of the receiver loops required to 
acquire and track the ranging signal components will involve an 
unavoidable subsystem design compromise. Wide tracking bandwidths 
are required to reproduce the large doppler dynamics caused by 
the spacecraft orbital motion on the high frequencies employed, yet 
narrow loop bandwidths are also desired to restrict the additive 
noise filtered along with the signal phase fluctuations. The dynamic 
range of the received signal level and the effects of corresponding 
SNR variations on the tracking performance must also be accounted 
for in the subsystem design to assure above- threshold operation and 
smooth departures from nominal operating conditions. 

The specification of the range and range rate data extraction 
and processing following the received signal reproduction also re- 
presents a significant subsystem design consideration. The range 
extractor must process the high-accuracy and ambiguity-resolving 
components of the range signal, so as to efficiently provide the 
range data without compromising tracking performance insofar as 
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possible. The range rate data extraction can be conveniently based 
on dcppler cycle -counting methods, which have proven their effective- 
ness in space tracking applications. The design parameters Involved 
in the measurement and arithmetic units of thc3e data extractors 
must be carefully chosen to minimize the bias and random error con- 
tributions to the range and range rate samples extracted. The out- 
put data rates and formats must, in turn, be compatible with the 

* 

pre-processing smoothing and orbital determination programs employed 
to track the spacecrafts in question. 
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3 . Ranging Signal Alternatives 

The ranging signal design consists of first identifying the 
ranging signal event and the doppler signal carrier respectively 
used to derive the independent range and range rate data , and 
then specifying the modulation principles employed to incorporate 
their characteristics into a unified ranging signal structure. 

The quantitative choice of modulation parameters and frequencies 
can then be based on the accuracy performance requirements of the 
application in question. There are two basic ranging signal de- 
sign approaches that are employed in range and range rate tracking 
systems: harmonic sidetone ranging ard pseudo-noise code ranging. 
These ranging signal alternatives are next discussed, along with 
their ranging system design implications and joint application, 
capabilities . 

3 . 1 Harmonic Sidetone Ranging Signals 

The use of continuous sinusoidal waveforms for ranging pur- 
poses is based on identifying the ranging event as the zero-crossing 
of the sinusoidal pulses. The period (T,) of the sinusoid then 
establishes the maximum unambiguous distance that can be resolved as 
T,j in two-way ranging, where (c) is the speed of light. The 
use of large periods (low frequencies) is incompatible with high- 
accuracy requirements, since range errors contributed by random 
noise effects are inversely proportional to the tone frequency (f,). 
For example, a random phase error (A4>) rad in the tone reproduced 
at the ground receiver corresponds to a range error Lr * 
in two-way ranging. The range error will thus increase with lower 
frequencies (larger periods) if the phase error in question is 


/c\ M> 

\ 2 / 2lTf ‘ 


- 11 - 


frequency- independent , which is the case for random noise error 
contributions . 

The design compromise between maximum ranging distance and 
high-accuracy performance is resolved by including other harmoni- 
cally related sinusoids of lower frequencies in the ranging signal. 
The highest frequency tone (major tone) is then used to provide 
small random error contributions to the range data extracted, 
while the remaining tones (minor tones) are used to resolve the 

t 

range measurement ambiguities caused by the major tone period. 

The ambiguity resolving tones are usually harmonically related, 
i.e., each tone frequency is an integer factor smaller than the 
preceding one of higher frequency. Every lower frequency tone is 
then used to resolve the ranging ambiguities of the preceding one 
of higher frequency, with the lowest frequency tore being capable 
of performing an unambiguous range measurement by itself. This 
last tone would yield an excessive randoL a rror contribution if 
used to provide the range data, but its role is simply to resolve 
the ranging ambiguities of the preceding tone so that its error 
tolerance is considerably relaxed. 

The use of only two tones is theoretically possible: a high- 
frequency tone (f h ) for accuracy resolution, and a low-frequency 
tone (f,) for ambiguity resolution. The practical problem lies in 
the tight phase error tolerance of the latter if expected to re- 
solve the ranging ambiguities of the former. Indeed, the maximum 
time error permitted for the low tone is a period (l/f h ) of the 
high tone, so that the low tone phase error magnitude is bounded 
by (2"f,) (1/f,,) - 2Tt(f, /fj rad. The frequency ratio (f h /f, ) has 
a large value for the high accuracy and long distance requirements 
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of spacecraft tracking, and the phase error bound would be exceeded 
in practice with just two ranging tones. The use of successive 
minor tones as explained above then relaxes the phase error bound, 
since the time error tolerance would then be the period of the next 
tone rather than that of the major tone. The frequency ratio and 
total number of minor tones employed is thus governed by the phase 
error performance in the reproduction of each tone at the ground 
receiver, so as to permit the successful resolution of the ranging 
ambiguities of the preceding tone. 

The mean-square phase error contribution from random noise 
effects is inversely proportional to the SNR existing in the track- 
ing loop reproducing the major tone. The presence of the minor 

tones thus limits the range accuracy performance by taking power 
away from the major tone loop SNR. This effect can be somewhat 
reduced by choosing modulation parameters that result in an optimum 
power sharing between the major and minor tones to achieve the best 
compromise in range accuracy vs. ambiguity resolution. However, 
the sequential (rather than simultaneous) transmission of the tones 
completely eliminates this power sharing compromise and shall be 
asiv’ied. The ambiguity resolution and fine ranging operations 
then occur sequentially, and only the major tone modulates the 
carrier during the fine range measurement. Of course, the power 
sharing between the tone (range signal) and carrier (range rate 
signal) remains as an accuracy performance compromise. 

The major tone frequency employed is selected from accuracy 
considerations as previously explained, with higher frequencies 
yielding smaller random measurement errors for given loop SNR con- 
ditions. However, the use of higher frequencies is accompanied by 
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larger doppler dynamics and oscillator noise effects, which re- 
quire wider acquisition and tracking bandwidths . These larger 
bandwidths then reduce the loop SNR and increase the random noise 
error contribution. An optimum loop design should thus jointly 
minimize the dynamic lag, oscillator noise and random noise phase 
errors in the reproduced major tone. Yet the loop optimization 
process must account for the different statistical nature of the 
phase error processes involved, so as to distinguish between bias 
and random errors in the range data extracted. An analogous dis- 
cussion applies to the minor tone loop design, except that the 
phase error tolerance is established from ambiguity resolution 
rather than range accuracy considerations, 

The lowest minor tone frequency employed is selected to un- 
ambiguously resolve the maximum range expected, and can result 
in rather low modulating frequencies. The modulation of these 
low frequencies on the signal carrier can result in signal 
components rather close to the carrier frequency, which may 
hinder the desired carrier extraction and ranging signal demodula- 
tion at the ground receiver. This situation can be avoided by 
exhibiting the lowest tone as a complemented bias on another fre- 
quency suitable for the signal generation and extraction operations. 
The use of SSB modulation on another higher frequency tone acting 
as a subcarrier has been used effectively in ranging system design 
under the terminology of doubly-complemented sidetones. This 
technique is applicable to any minor tone whose frequency is too 
low for a reliable system processing without double complementa- 
tion. 
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3 . 2 Pseudo-Noise Code Ranging Signals 

The ranging event in pseudo -random signal consists of the wave- 
form structure itself, and the range measurement is based on the 
time shift required to match the transmitted and received waveforms. 
The notion of matching requires a characteristic property or para- 
meter associated with the waveform that provides a reliable recog- 
nition of zero delay, plus easily discriminates when time shifted 
replicas are present. The use of the waveform auto-correlation 
function for time delay discrimination represents a natural choice 
since its magnitude peaks at the origin. An arbitrary analog pseudo- 
random waveform with a peaked auto-correlation function is perfectly 
valid, but it may introduce excessive storage and processing require- 
ments, since the waveform duration (repetition period) is selected 
to provide unambiguous ranging to the spacecraft. The use of a 
coded waveform relaxes this consideration, since the number of digits 
and not the period in seconds governs the storage and processing 
complexity. 

The use of pseudo-noise codes is thus motivated as an alternative 
ranging signal. Their autocorrelation function at integer shifts 
of a symbol duration is given by 


R (t) = 


1 

1 

L 


for t = k P 
for t * k P 


where (k) is an integer, (p) is the code period, and L = (2" - 1) 
is the code length or number of code digits in a period, Such 
auto-correlation function can easily discriminate the matched and 
mismatched states between transmitted and received codes. The 
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code period (P) is selected to provide the unambiguous ranging, 
while the symbol duration (P/L) determines the code bandwidth 
occupancy. The use of short symbols is desired since the range 
errors contributed by random noise effects are inversely propor- 
tional to the symbol rate (L/P) . The latter plays a role analo- 
gous to the tone frequency in error control, but without any 
ambiguity resolution limitations. Indeed, the use of high symbol 
rates (short symbols) and long code periods are simultaneously 
possible by having a large code length (L) . This represents an 
extra degree of freedom relative to sidetone ranging, where the 
tone frequency and the ranging event period were essentially the 
same parameter for design purposes. 

The code length (L) governs the duration of its acquisition 
or matching process. The maximum number of (L) trials required 
can be reduced by replacing the code by a boolean combination of 
subcodes of smaller relatively prime lengths (Lj). The resultant 
code length is given bv the sum L = X Li , while the resultant code 
period is given by the product P = n Pi. The original code period 
and unambiguous range capability can thus be maintained with a 
substantial reduction in the maximum number of test trials involved 
by acquiring the subcodes separately. However, each subcode 
matching process will not exhibit the same detection immunity present 
with the original code autocorrelation function, and therein lies 
a compromise of acquisition time vs. code detectability in the 
presence of noise. 

The repetition of the code matching process for consecutive 
range measurements can be avoided by tracking the code after acqui- 
sition. This may b< implemented with delay-locked loops which track 


- 16 - 


a demodulated code, or with code-clock double loops that also ex- 
tract a coherent clock component. Both these processors represent 
tracking cross-correlator realizations analogous to the conven- 
tional phase-locked loops employed to track the sinusoidal tones. 
However, the simpler processing and reliable experience available 
with the conventional loops represents an advantage of tone ranging 
over code ranging from acquisition and tracking considerations. 

The long-range ambiguity resolution advantages of code ranging 
can be maintained without sacrificing the desired tone ranging fea- 
tures by having the code replace the lower frequency tones. This 
hybrid ranging signal then derives the fine range measurement from 
the major tone, and resolves range ambiguities with some interme- 
diate minor tones plus the code. One of these minor tones can also 
be used as a code subcarrier to derive a clock component. The 
reduction in the number of minor tones employed when compared to the 
pure sidetone case is particularly desirable in our application, 
where the large major tone frequency required from accuracy consid- 
erations would induce a large number of such minor tones to effec- 
tively resolve all ambiguities. 


4. Signal Modulation and Processing 

The use of harmonic sidetone and/or pseudo-noise code as 
ranging signals represents an effective approach towards an inde- 
pendent range and range* rate system design. A hybrid ranging 
signal composition is most advantageous to exploit both the tone 
and code signal features. The sequential application of the fine 
ranging and ambiguity resolving operations is also desirable to 
minimize the power sharing compromise between the range and range 
rate signal components. Under these conditions, the ranging signal 
structure during the fine ranging operation consists of the laser 
carrier modulated by a single ranging tone whose frequency is 
large enough to meet the range accuracy requirements . The doppler 
effect in the received carrier then provides the range rate data, 
and the phase offsets in the demodulated tone provide the range 
data independently. 

The modulation of the ranging tone on the carrier presents 
various choices to the system designer. The use of amplitude 
modulation or narrowband phase modulation results in a double 
sideband spectrum around the carrier frequency, without any base- 
band signal bandwidth expansion. The use of a single sideband 
spectrum would further reduce the channel bandwidth, but the 
demodulation process may introduce a severe degradation in the 
range accuracy performance as illustrated below. The use of wide- 
band phase modulation is not effective since it expands channel 
bandwidth while setting a limitation on carrier power allocation. 
The use of coupling modulation represents another alternative in 
laser ranging systems, and its characteristics are also discussed 
below. 
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The use of AM techniques can be considered with the provision 
that a carrier component is mi intained to permit the independent 
range rate data extraction at the receiver. The AM signal without 
carrier suppression can be formulated as 

v AM (t) = A [l + a m(t)j cos (u> c t + e c ) (A. la) 

= A cos (w c t + 0 C )± ^ cos £(u> c ± ww)t + 0 C ± 6 m j 

for m(t) = sinCui^t + 0.) (A. lb) 

The modulation factor (M.F.) is defined as (a) times the peak 
value of the modulation m(t) , which reduces to M.F. * a in our 
sinusoidal case. In turn, the peak factor (P.F.) is defined as the 
ratio of the peak to rms values of the modulation m(t), which re- 
duces to P.F. = \fl in our case. The average power in the AM 
signal can then be evaluated in terms of these factors to yield 

v»(t) = | ! i + a 2 m 2 (t) = f 2 [i + (ftf :) 2 <4 2a) 

= | (1 + | ) for m(t) = sin (u) m t + 0 m ) (A. 2b) 

where the derivation assumes modulation frequencies well below the 
carrier frequency, plus zero d-c value in the modulation. Under 
these conditions, the total average power decomposes into a carrier 
term ^ plus a sideband term ^ f -1 ) ‘ The ratio (M.F. /P.F.) 

then governs the power allocation between the sideband (range) and 
the carrier (range rate) components in the ranging signal. Such 
ratio reduces to a 2 /2 in our sinusoidal modulation case, so that 
the M.F. = a becomes the power sharing design parameter. 
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The use of envelope demodulation for the modulation extraction 
at the receiver requires M.F. < 1, so as to avoid carrier phase 
reversal and permit the modulation recovery. The conditions 
M.F. < 1 and P.F. > 1 then illustrate that no more than 50% of the 
total power can be allocated to the sidebands in an AM signal with 
envelope detection. In the sinusoidal modulation case, the con- 
ditions become M.F. = a < 1 and P.F. = JY, so that no more than 

ii 

33% of the total power is available to the sideband (range) com- 
ponents. These constraints on the power allocation between the 
sidebands (range) and carrier (range rate) components are a direct 
consequence of the envelope demodulation assumed. The interest 
in such receiver processing is further diminished by noting that 
the carrier extraction is anyway required for range rate purposes 
regardless of the tone demodulation process. The carrier extraction 
via phase-locked techniques ideally yields a coherent carrier 
replica at the receiver, which can then be used for product demodu- 
lation purposes. This approach does not impose any a-priori condi- 
tions on the modulation factor, hence relaxing the power sharing 
limitation inherent in envelope demodulation. 

The product demodulation of the AM signal with a coherent 
carrier yields the modulation at the receiver. In the sinusoidal 
modulation case, the demodulated tone may be jointly tracked and 
filtered with a phase-locked loop to permit continuous range estima- 
tion. The modulation parameter (a) is then selected to provide 
enough sideband SNR for the tone phase-locked extraction and accurate 
range measurement, while maintaining enough carrier SNR for the 
carrier phase-locked extraction and accurate range rate measurement. 
The quantitative analysis should account for the 3db improvement 
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in the tone SNR due to the coherent sideband folding, as well as 
the degradation introduced by the lack of ideal coherence in the 
reproduced carrier replica. A carrier locking error (A<J>) causes 
an amplitude attenuation by cos(A<j>) in the demodulated tone, with 
a corresponding SNR degradation by a cos 2 (A<f>) factor. 

4. 2 Angle Modulation Systems 

The constant envelope features of PM or FM techniques have 
motivated their use in existing range and range rate tracking 
systems. The AM signaling approach may require a low-level opera- 
tion to avoid nonlinear distortion or saturation effects, whereas 
angle modulation permits efficient nonlinear amplification provided 
AM/PM conversions are controlled. The distinction between PM and 
FM systems is particularly absent in the narrowband sinusoidal 
modulation case, using coherent product demodulation at the receiver 
for the tone extraction. 

The use of a low modulation index (3) in the PM signal estab- 
lishes a correspondence with the AM case when product demodulation 
is employed. The PM signal may be formulated: 

v fM (t) = A cos ju) c t + 0 C + 3m(t)”J (4.3) 

* A cos (o) c t + 6 C ) - RAm(t) cos (<u c t + 6 C ) for B << 1 (4.4a 

(o) c ± cO t + e c ±0. 

for 6 << 1, m(t) - sin (w w t +6.) (4.4b) 

1 

so that a comparison of (4.4a) and (4.4b) with the corresponding 
AM signals expressions in (4.1a) and (4.1b) shows that the only 
difference is a quadrature carrier in the PM case vs a cophasal 
carrier in the AM case. Such distinction does not affect the power 


(w c t + 0 C ) + y sin 
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content, so that the results of (4.2a) and (4,2b) are also applica- 
ble with the paramters (B) and (a) interchanged. The disadvantage 
of low- index PM signaling is then that the power allocated to the 
sideband (range) components is a small fraction of the total power 
available . 

The use of arbitrary index values does not have a general 
compact formulation, but for sinusoidal modulation it can pre- 
serve the narrowband spectrum (carrier plus one sideband pair) with- 
out introducing extreme power sharing constraints. The PM signal can 
be decomposed as a carrier plus an infinity of sideband pairs as 


v rM (t) = A cos |"u) c t + 9 C + B sin (u)*t + 9 m )j 

oo 

= A J n (B) cos ± na, » t ) + 0 c ± n0 » J 


(4.5a) 

(4.5b) 


n =-oo 

where image sidebands have the same magnitude since the Bessel 
coefficients satisfy the symmetric relation J_ n (B) * (-l) n J n (B). 

The total average signal power is always given by A 2 / 2 regardless 

00 

of the index value, since ^ ^ (B) ■ 1 for all B. Each individual 

n= - oo 

sideband component power is down by a factor J n 2 (B) from the total 
power, so that the modulation index (B) governs the power sharing 
distribution between the carrier and all the sidebands. The use 
of 0 < B <1.5 further assures that only the carrier (n = 0) and 
the first sideband pair have a significant power allocation as 
illustrated by Figure 1 and Table 1. This case results In a narrow- 
band spectrum without the B << 1 constraint as a consequence of the 
Bessel functions involved. The use of larger modulation indeces 
would waste power in higher order sidebands, so that wideband 
angle modulation is not considered. 
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TABLE 1 

PM Carrier vs Sideband Power Sharing 
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0.1 
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98.01 
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0.01 

0.3 
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4.39 
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0.4 

92.24 

7.68 

0.08 

0.5 

88.08 

11.74 

0.18 

0.6 

83.17 

16.44 

0.39 

0.7 

77.65 

21.65 

0.70 

0.8 

71.62 

27.20 

1.18 

0.9 

65.21 

32.95 

1.84 

1.0 

58.55 

38.73 

2.72 

1.1 

51.78 

44.35 

3.87 

1.2 

45.04 

49.66 

5.30 

1.3 

38.46 

54.49 

7.05 

1.4 

32.14 

58.73 

9.13 

1.5 

26.19 

62.25 

11.56 

1.6 

20.73 

64.96 

14.31 

1.7 

15.84 

66.77 

17.39 

1.8 

11.56 

67.63 

20.81 

1.9 

7.94 

67.56 

24.50 

2.0 

5.01 

66.52 

28.47 
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The tone demodulation is again performed via phase- locked 
carrier extraction followed by coherent product demodulation and 
phase-locked tone filtering. The only distinction with the pro- 
duct demodulation operation of AM signals consists of a quadrature 
phase shift on the reproduced carrier to derive a coherent carrier 
replica. The modulation index (3) is selected to provide enough 
carrier and tone SNR for locking purposes, as well as accurate 
range and range rate data. Again, the 3db SNR processing gain, 
and the cos 2 (A$) tone SNR degradation in the presence of carrier 
locking errors (A<t>) have to be accounted for in the analysis. The 
use of narrowband PM can thus match the signal bandwidth, power 
sharing and receiver detection characteristics of AM with coherent 
product demodulation, plus provide a constant envelope signal for 
effective processing. 

A. 3 Single Sideband Systems 

The spectral occupancy of the ranging signal can be minimized 
by using SSB signals. Such signals can be obtained by either 
filtering the undesired sideband from a DSB spectrum, or by the 
phase shift method of SSB signal generation. The latter is based 
on the recognition that a SSB signal is a joint AM-PM signal, which 
can be decomposed as the sum of two AM-DSB signals on quadrature 
carriers. The two quadrature modulations are related via a Hilbert 
Transform, which simplifies to a simple 90° phase shift in the 
sinusoidal modulation case. The net effect of the dual modulation 
is the cancellation of the undeaired sideband component. 

The tone demodulation at the receiver proceeds as before 
through coherent product demodulation. There is no 3 db SNR 


processing gain since no sideband folding occurs, but this is 
exactl ’ compensated by the extra amount of power that can be put 
in thu >in»le sideband component since the other sideband is not 
transmitted. However, the presence of a carrier locking error 
will produce rather different and harmful effects when compared 
to the DSB case. The product demodulation will now transfer the 

carrier phase error (A<|>) to the demodulated tone as an additive 

*• 

phase shift, rather than as the amplitude attenuation effect ob- 
tained in the DSB case. 

A comparison of both effects is illustrated in Figure 2, The 
SSB case introduces a time error of At = in the reproduced tone 

U) m 

zero crossing, with a corresponding range error of Ar = (A<p) . 

2w m 

Notice that these errors are essentially magnified by the carrier/ 
tone frequency ratio since the phase error (A<|>) pertains to the 
carrier reproduction, yet it is divided by the tone frequency (not 
the carrier frequency) to produce the time and range errors . The 
fact that both carrier and tone locking errors appear as additive 
range errors via the same conversion factor is a considerable 
degradation relative to the DSB case. The latter limited the 
carrier effect to a decrease in the tone SNR. by cos 2 (A<J>), which 
only increases the rms error contribution of thermal noise alone 
by a sec(A<p ) factor, but does not introduce any additive magnified 
error effects. The use of SSB signals is not employed in indepen- 
dent range and range rate tracking systems on this basis. 

4 . 4 Coupling Modulation Systems 

The previous approaches have always identified the range tone 
as the modulating function on the carrier. The modulation process 
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(1) DSB : cos (A4>) sin (oj m t + 0m) 
SSB : sin (<i) m t + 6 m + A<J>) 


(2) DSB: 
SSB: 


sin (o) m t + 0„ + A0m ) 
sin (oj m t + 0™ + A0 * + A<J>) 


( 3 ) 


DSB: rms range error due to additive noise increases by a sec (A<|>) multiplicative factor 
SSB: total range error increases by a A<t> additive term. 


Figure 2: Carrier Locking Error Effects in DSB vs. 


SSB Systems 
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produced a carrier plus first-order sidebands as the desired spectrum, 
with the tone frequency appearing as the differential frequency be- 
tween carrier and sidebands. Such spectrum could also have been 
obtained by having the modulating function be another periodic 
function containing the tone frequency as the fundamental or a har- 
monic of it. For example, a rectangular pulse train with its funda- 
mental as the tone produces a PM spectrum that can be decomposed in 
the form of (4.5b), except that the coefficients are now given by: 

C (6) = sine ^ (3 - n) instead of the J (3) Bessel coefficients, 
n d Tii 2 n 

The use of coupling modulation is representative of such approach. 
It essentially consists of an AM system, where the modulating function 
consists of a sinusoidal functional of the tone frequency. The 
modulated waveform can be decomposed as a carrier plus sidebands as 
follows : 


v cm (t) = ^ sin ^Y 0 + Y J sin (w^t + 6 m )J cos Cu) c t+e c ) (4.6a) 


= A 0 (sin Yo) J Q (Yi) cos (u) c t + 0 C ) 


00 

A o ( s in Yo) 2J n (Yi) cos ^(u> c : nu tn ) t + 9 C t ne*, 1 

, even 

r oo 

Ao (cos Yo) ^ 2J m (Yi) cos. j^(oo c t nmj t + 9 C - n9 m J(4.6b) 


The power sharing distribution is governed by the two modulation 
parameters (Yo) and (Yi). The case of Yo = 0 yields a suppressed 
carrier, which is undersired in an independent range and rate rate 


system, plus prevents coherent product demodulation even if no range 
rate data is to be extracted. The case of Y 0 = ~ suppresses all 

odd sidebands, so that the logical power sharing compromise involves 
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( Y i ) for the carrier vs. 2j| (Yi) for the n - 2 sidebands. This 

case yields an inferior power sharing option when compared to the 

2 2 

J Q (B) vs. 2J^ (6) compromise existing in PM systems. The minimum 
bandwidth spectral allocation of a carrier plus first order (n = 1) 

sidebands results in a power sharing compromise of J Q ( Y 4 ) sin 2 Y 0 vs. 

2 

2J^ ( Y x ) cos Y 0 . which again represents an inferior option relative 
to PM systems. In any case, it should be noted that coupling modula- 

t 

tion does not yield a constant total power waveform, since the rela- 

00 

tion (Yi) +2 ^ (Yj) = 1, cannot be applied; i.e., the total 
average signal power is given by 


- | { (sin 2 Y 0 ) 
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(4.7) 


4. 5 Laser Modulation Specifications 

The previous discussions on the ranging signal structure and the 
power sharing compromises involved should also account for the 
allowable parameter values. The laser specifications limit the 
modulation frequency to 5 MHz for AM and 20 MHz for PM, and the 
channel bandwidth to 300 MHz. The modulation index is also restricted 
to 1 rad maximum in the PM case. The lower modulation frequency 
allowed for AM will restrict its random range error accuracy perfor- 
mance , since the latter has rms values inversely proportional to the 
tone frequency. However, a higher modulation frequency will imply 
larger doppler dynamics to be tracked by the tone loop, hence wider 
loop bandwidths and increased phase noise locking errors due to 
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additive thermal noise. The error analysis should thus carefully 
account for both locking and accuracy limitations when selecting the 
tone frequency. A quantitative evaluation is required to establish 
which of these two effects plays a more predominate role. 

The channel bandwidth limiation is a secondary effect since 
narrowband conditions are present in AM and are desired in PM for 
effective power sharing. The modulation index restriction in PM is 
a rather important specification since it limits tone power allocation 
(and range accuracy performance) as previously illustrated in Table 1. 
The maximum index value of 1 radian would assign 58.557o of the 
available total power to the carrier, 38.73% to the tone sidebands, 
and 2.72% to the higher order sidebands. An increase of the allow- 
able index value to 1.4 radian would essentially reverse the alloca- 
tion to 32.14% to the carrier, 58.73% to the tone sidebands, and 
7.05% to the higher order sidebands. Such an extension of the 
modulation index range is noted to result in a significant improve- 
ment in the power sharing options available in PM systems . The 
quantitative analysis shall include such extended range to project 
error performance capabilities of the laser ranging system under 
consideration . 
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5 . Phase-Locked Loop Tracking Performance 

The ranging signal spectrum essentially consists of a carrier 
component plus a range tone sideband pair during the fine ranging 
operation. The receiver then phase locks to the carrier component 
and uses it to coherently product demodulate the range tone. The 
latter is also tracked with a phase locked loop and used to pro- 
vide range data, while the extracted carrier provides independent 

» 

range rate data. It should be emphasized that even if no indepen- 
dent range rate data is to be derived, the carrier component must 
have sufficient power to permit carrier phase lock in the presence 
of additive noise, so as to assure an effective coherent demodula- 
tion of the range tone. Hence, the carrier power allocation is 
established from carrier lock and range rate accuracy considerations 
(whichever predominates), while the tone power allocation is esta- 
blished from tone lock and range accuracy considerations (whichever 
predominates) . 

The phase locked loop designs must account for both the acqui- 
sition and tracking operations. The acquisition design must assure 
that the loop can acquire the maximum doppler plus drift conditions 
expected in the presence of additive noise with a low probability 
of failure in a moderate time. The design is based on existing 
data on the probability of acquisition as a function of the loop 
SNR, loop parameters and initial frequency offset to be captured. 

The use of manual or automatic acquisition aids such as VCO sweeps 
can be exploited to considerably ease the acquisition process and 
relax its loop design compromises. Once the signal acquisition 
has been successfully accomplished, the loops behave as equivalent 
linear time invariant filters on the input signal phase fluctua- 
tions plus additive or oscillator noise effects. The tracking 
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design is then based on maintaining the locked state with a small 
phase error. The loop design involves a basic compromise: wideband 
loops are effective in tracking doppler dynamic and oscillator jitter 
phase inputs, yet narrowband loops are desired for additive phase 
noise rejection. The phase error contributions of these super- 
imposed effects should maintain the locked state and validate the 
linear filter model. 

The baseband block diagram and linear phase transfer model 
applicable to phase-locked loops is shown in Figures 3(a) and 
3(b). The input signal s(t) exhibits a time-varying phase 4>(t) 
assumed to induce all oscillator jitter (including the VCO 
effects) besides the doppler phase dynamics to be reproduced. 

The VCC reference signal r(t) produces an estimate of such input 
phase in its output phase <j> (t) , and the phase detector develops 
an error signal u(t) which is a (2tt) periodic function f(4>) of 
the phase estimation error 4 > ( t ) = <f\ (t) - <J> o (t) , The additive 
noise n(t) accompanying the input signal is decomposed as cophasal 
and quadrature components about the nominal signal phase, and it 
contributes an additive noise component n'(t) at the phase detector 
output. The loop filter F(s) produces a VCO control signal v(t) 
which ideally causes a proportionate change in the VCO offset fre- 
quency 4> 0 (t) from its free-running value. 

The small phase error conditions prevailing in the locked state 
permit the linearization of the phase detector characteristic 
f($) ss <J> to obtain the linear phase transfer model. The output 
phase estimate $ G (s) and phase error ^(s) transforms are then 
related to the signal phase 4> . (s ) and equivalent phase noise 
$ (s) input transforms through the transfer functions 
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Figure 3: Phase-Locked Loop Block Diagram and Linear Models 
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$ . (s) + <t> 



- 

KF - ^ 4> (s) + <j> (s) (5.1) 

s + KF(s)L J 


D (s) « [\ - H ( s )] $ , (s) - H(s) *„(«) 

= - 4>,(s) - « (s) (5.2) 

s + KF (s) s + KF(s) " 

The closed-loop transfer function H(s) represents the equivalent 
filter estimating the signal phase fluctuations in the presence 
of additive noise. The loop design consists of specifying the 
loop gain K and loop filter F(s) to yield an effective phase esti- 
mation. The first error term in (5.2) consists of the inability 
of the loop to ideally reproduce the signal phase fluctuations, 
and it includes the dynamic lag and oscillator noise errors. The 
second error term in (5.2) consists of the extra phase error 
contributed by the additive noise passed by the equivalent loop 
filter. The loop design compromise is illustrated in this formula- 
tion: a wideband H(s) reproduces the high frequency content in 
<I> j ( s ) , but a narrowband H(s) rejects the high frequencies in 
$n(s) • 

The steady- state dynamic lag error when tracking polynomial 

step excitations of the form <J> f (t) = — M t" for t > 0 can be 

1 n! n 

evaluated from the final value theorem as 

<t> = lim s [l - H(s)l <D (s) » lim S ' (5.3) 

Such input phases may be obtained from a power series expansion of 
the orbital dynamics and corresponding doppler dynamics; so that 
Mo represents a constant phase offset (constant range) , H, is a 
constant doppler shift (constant velocity) , M 2 is a constant doppler 
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so that (m + 1) - th order phase dynamics could be tracked with 
finite steady-state error (and m-th order with zero error) if (m) 
ideal integrators are present in the loop filter. The non-ideal 
integrator case can only track 2nd-order phase dynamics (constant 
doppler shift) with finite steady-state error (and constant phase 
offsets with zero error) . The distinction is rather important 
since it invalidates the use of a steady-state error analysis based 
on ideal integrators, when the latter are recognized as absent in 
practical loop filter realizations. 

The use of a steady-state error analysis is questionable by 
itself. The power series expansion is valid about a point in the 
orbital motion, so that letting t — *oo may be incongrous since the 
orbital dynamics char 0 t with time. For example, the non-ideal 
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integrator result indicates that a constant doppler rate cannot 
be tracked forever without a time-increasing phase error, so that 
the locked state and the loop linear model would be lost. Yet 
this does not exclude that such doppler rate cannot be tracked 
for a finite amount of time within the locked state, which may 
suffice if the tracking time in question is larger than the 

orbital time that validates the constant doppler rate assumption. 

* 

As time elapses, the doppler rate does not remain constant, and the 
dynamic tracking error should be reformulated. Notice also that 
tracking a peak doppler rate magnitude for a short time is not 
necessarily more difficult than tracking a moderate magnitude 
for a longer time. 

The use of an rms error analysis for the doppler dynamic 
phase inputs under consideration should be avoided. The rms evalua- 
tion consists of the (magnitude-square) time -averaging of the 
transient error response to these phase input excitations. Such 
error averaging includes transient overshoots, damped oscillations 
and steady-state error sections. An infinite steady-state error 
would always yield an infinite rms error, so finite error magni- 
tudes would be limited to ideal integrator designs, or non-ideal 
integrators with constant phase offset or doppler shift inputs. 

Also, a finite rms error derived from the linear model analysis 
cannot assure that the error peaks present in the transient over- 
shoot and ringing phenomena are within the locking bound (which 
would invalidate the linear model and the error analysis). Finally, 
the transient time-averaging would refer to orbital time, yet the 
rms error contribution of additive noise represents a stochastic 
ensemble averaging applicable to only one orbital point (since the 
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path loss fluctuation with orbital motion would chdnge the SNR con- 
ditions) so that the Joint accounting of rms dynamic plus noise 
effects would suffer interpretation. 

The rms additive noise error follows from (5.2) in terms of 
the power density spectrum S (oj) of the equivalent phase noise input 


(add) 


V/„ v- : 




H(jo) 


) 


2 2 
2tt rad 


(5.6 


The power density spectrum S n * (w) of the phase detector baseband 
noise n' (t) may be evaluated in terms of that of the input noise 
as a function of the VCO output phase fluctuations. If the closed- 
loop equivalent filter is narrowband relative to the loop input 
passband, the VCO phase effect will approximate an orthogonal 
transformation that preserves the baseband noise spectra. Under 
such conditions, the power density spectrum of n' (t) is similar to 
those of n c (t) and n (t) ; so that a two-sided N q watts/Hz white 
spectrum may be assumed when the noise input n(t) has a one-sided 
N q watts/Hz white spectrum. The rms additive noise error can then 
be expressed in terms of the two-sided loop noise bandwidth 

(B n in Hz) and the corresponding loop SNR via 

/oo 

|h(jw)| 2 Hz (5.7) 
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The rms oscillator noise error also follows from (5.2) by 
using the power density spectrum (u) of the equivalent oscillator 

i 

phase noise at the loop input to ootain the relation 
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Such error contribution cannot be given a more compact formulation 
without further a-priori information, since the oscillator noise 
processes superimposed at the loop input depend on both the signal 
sources and subsystem processing employed in a given application. 

The selection of loop gain K and loop filter F(s) control 
these error magnitudes only through their effects on the equivalent 
transfer functions H(s) and 1 - H(s). The latter specify loop order 
as the degree of their denominator characteristic polynomial in the 
(s) frequency domain. The next subsections shall consider loop 
filter realizations that result in second or third order loops. 

Their dynamic lag and additive noise error contributions shall be 
analyzed towards an optimum loop design, with the oscillator noise 
effects yet to be accounted for. Once these tracking errors are 
formulated, their corresponding range and range rate error contri- 
butions can be evaluated in the sections that follow. 

The set of loop design parameters thus consists of the loop 
gain and the loop filter time constants or transfer functions. 

However, the following two results are applicable to all loop orders 
and filter realizations that satisfy F(0) - 1 (non-ideal integrators): 

(a) The steady-state dynamic lag error when tracking 

a constant frequency offset of ft rad/sec is essentially 

n 

given by rad in the locked state (see (5.4)); i.e., 

I\ 

f " * ( k ) ® it for errors. 
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(b) The rms additive noise error in the locked state 
is given by (5.7) in the presence of wideband additive 
white noise. 

The importance of these results is that different loops will exhibit 
the same steady-state dynamic lag error (to a constant frequency 
offset) and the same rms additive noise locking error provided they 
have the same loop gain (K) and looo noise bandwidth (B n ), regard- 
less of the other loop parameter values and the loop orders in 
question. 

5. 1 Second-Order Loops 

The class of 2nd-order loops under consideration is character- 
ized by the presence of a passive or active lag network as their 
loop filters. The loop filter realizations are illustrated in 
Figure 4, along with tueir corresponding transfer functions. The 
active case provides an increase in the realizable pole time con- 
stant with the operational amplifier gain as the magnification 
factor. Both the passive and active realizations have been employed 
in existing tracking systems. The existence cf non-ideal integra- 
tors is emphasized in both cases, with the ideal integrator being 
closer approximated in the active case due to the inclusion of the 
amplifier gain in the time constant. The loop filter and closed- 
loop transfer functions can in both cases be formulated as 

F(s) - 1 + - T ^ , T < (5.10) 

1 + T 2 S 1 2 

and 

H(s) . i j T > 8 

1 + Ci + h s + li s ! (5.11) 

K K 
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The set of loop design parameters may then be taken as the 
trio (K, ij , x 2 ). The alternate set (K, u r , O is more often 
employed based on servo-loop analogies, where (oj r ) is the reso- 
nance frequency in rad/ sec and (?) is the damping factor specified 
by 

-f" . 5 ■| (T . + | ) “r (512) 

&J 

The equivalent (two-sided) loop noise bandwidth (B n in Hz) can 

then be evaluated in terms of either set to yield 

) 4c<J r 

. 1 + - T^ + (ifj Ur (5.13a) 

Vk 

( — + ~ | = 1 + 4c * u r for Kt, » 1. K » to (5.13b) 

2 \t, t 2 ) 4c r 

The variation of the loop noise bandwidth (B n ) with the damping 
factor (s) is quite small for 0.5 ^ ^ 1, which represents con- 

ventional 2nd-order loop designs; e.g., B n /w r = (1, 3/2 J~2, 5/4) 
for £ = (1/2, 1). It should be emphasized that the results 

above have the units of B n in Hz vs. w r in rad/sec; e.g., 

B p (Hz) = w r (rad/sec) for C = 1/2, Notice that the set (K, B n , O 
represents still another option for the design parameters in a 2nd- 
order loop. This last set is preferable to the others because it 
features the loop gain and loop bandwidth as the design parameters 
which was previously motivated. An equivalent set that maintains 
this feature consists of (K, B n , a), where the last parameter 
satisfies 0 < a <2, and is suggested by the formulation of (5.13b) 
as follows: 
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— a a B p , — « (2 - a) B n 

T l T 2 


or 


a a 


x , Kti 1 + 4; 2 

t 2 


for K >> a B. 


(5.14) 


(5.15) 


The set (K, B n , a) could have been directly obtained from 
the closed- loop transfer function of (5.11). The interest in 
exhibiting the loop noise bandwidth as ar. independent design para- 
meter suggests the transformation 


1 

a B n 


II = _i 

K b B n 2 


(5.16) 


which can be introduced in the expressions of (5.13b) to derive 
the following functional dependence for Kt x >> 1: 


b as a(2 - a) , K >> a B n 


(5.17) 


Notice that the approach consists of first formulating the time 
constants (i lf — ) of H(s) as inverse functions of the hand- 
width but including an exponent on the latter to match that of the 
s-polynomial term. The set (K, t , t 2 ) is thus transformed into 
the set (K, B p , a, b) , and then the dependence b(a) is acknowledged 
to obtain the set (K, B n , a) with three independent design para- 
meters. Notice that in reality (b) is dependent on all three of 
these independent parameters as obtained from (5.13a), but the 
use of the approximation Kt j >> 1 and its result (5.13b) yields 
the simplified one-parameter dependence of (5.17). It is evident 
that this approach can be generalized to higher-order loop transfer 
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functions to free the loop noise bandwidth as an independent design 
parameter, and therein lies the motivation for this set. 

It should be noted that the approximate results obtained under 
the condition Kt x >> 1 become exact when ideal integrators are 
assumed in the loop filter. The deletion of the (1) constant and 
the (s/K) term in the denominator of (5.10) and (5.11) respectively 

t 

correspond to this assumption, from which the expression of (5.13b) 
then follows as an exact result. Hence, the loop noise bandwidth 
and rms additive noise error are essentially similar in the ideal 
vs. non-ideal integrator cases. However, their steady-state dynamic 
lag errors are considerably different in nature as illustrated 
below. The error transfer function governing such performance in 
the non- ideal integrator case is given by 


1 - H(s) = (5.18) 

1+ (t + -)s + — s 1 2 
1 k' K 

while both (s/K) terms in the numerator and denominator are absent 
for the ideal integrator case. The (1/K) term may, indeed, be 
neglected in the denominator of (5.18) for Ktj >> 1, but not in its 
numerator where it appears alone and may play a significant role in 
the steady-state error behavior. 

The steady-state dynamic lag error of a 2nd-order loop with 
non-ideal integrators is ft/K rad for a constant doppler shift 
(ft rad/sec) but infinite for a constant doppler rate (ft rad/sec 2 ) 
or higher dynamic inputs. However, the steady-state error behavior 


- 44 - 


I 



can still be studied for any input by performing a long division of 
(5.18) to derive the pertinent error coefficients. The expansion 
and truncation process yields the following results for the two 
inputs under consideration: 

♦„ * | (5.19) 

and 



The doppler rate error expression is noted to include a 
linearly time increasing term which causes the infinite steady- 
state error result. Such term is directly contributed by the 
(s/K) term in the numerator of (5.18), and would be absent if 
ideal integrators were assumed. Hence, the ideal integrator 
case could track a constant doppler rate forever with a finite 
constant steady-state error, but the non-ideal integrator case 
would result in the unlocked state as the tracking time increases. 
The conditions of high loop gains and short-lasting doppler rates 
relative to the tracking times are required if the non-ideal 
integrator error of (5.20) is to be approximated by its constant 
term only. In any case, the constant term is minimized to the 
value ft/B n 2 by the design e = 1/ /2, a * 1 , B n (Hz) * w, (rad/sec), 
which shall be assumed for the 2nd-order loops in question. 
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5 . 2 Third-Order Loops 

The class of 3rd-order loops under consideration is essentially 
characterized by a cascade of two passive or active lag networks 
as their loop filters. These loop filter realizations are illustra- 
ted in Figure 5, along with their corresponding transfer functions. 
Again the active networks are shown to provide an extension of the 
realizable pole time constants. The all-active realization has not 
been employed in existing tracking systems to the best of our 
knowledge. It should be noted that all the cases under considera- 
tion yield two real-axis poles in the s-plane, yet the formulation 
that follows shall be generalized to include the possibility of 
complex conjugate roots resulting from a different loop filter 
realization. The loop filter and closed-loop transfer functions 
can then be formulated as 


F(s) 


1 -t- T,S + (T,S)* Tj < t4 
1 + T, S + (t 4 S) 2 


and 


H(s) 


1 + T,S + ( T j S ) 2 

1 + (t, + i )S - (Tj + Ii)s ! + S J 
K K. K 


(5.21) 


(5.22) 


The set of loop design parameters may be taken as (K, T lf r 2 , 
t 3 , r l4 ) / but such set would reflect the loop filter rather than the 
closed-loop performance. The interest is to choose a set that 
includes the pair (K, B n ) as previously indicated. The loop noise 
bandwidth may be evaluated using (5.7) to yield 



(5.23a) 
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Figure 5: 
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Third Order Loop Realizations 





(5.23b) 



where the approximate result of (5.23b) becomes exact if ideal 
integrators are assumed in the loop filter; i.e., the deletion 
of the (1) and (t 2 s) terms in the (5.21) denominator, and of the 
(1/K) and (t 2 s 2 /K) terms in the (5.22) denominator, yields (5.23b) 
as an exact result from (5.7). 

The parameter transformation of interest is motivated by the 

closed- loop transfer function of (5.22) along with Kt j >> 1 and 

2 

Kt 3 >> t 2 as follows: 






(5.24) 


which replaces the original set by the set (K, B n , t 2 , a, b, c) . 

The non- ideal integrator parameter t 2 has a negligible effect on 
the closed- loop transfer function and loop noise bandwidth expres- 
sion, but becomes relevant in the error transfer function 1 - H(s) 
and the steady-state dynamic error analysis. The trio (a, b, c) 
contains only two independent parameters so as to maintain the 
same degree of freedom as in the original set. The dependent 
parameter is actually a complicated expression of all the parameters 
in the new set, but the approximating conditions in (5.32b) establish 
,a two-parameter dependence c(a,b) as the positive real solutions of 
the quadratic equation 

c’ + |" |1 - (a + 2) 1 c - 2ab 2 « 0 (5.25) 
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with an equal sign holding if ideal integrators are assumed. The 
situation is analogous to the pair (a,b) of 2nd-order loops, where 
the relation b « 2 - a established their dependence. The equation 
(5.25) can be noted to be quadratic for either of the three para- 
meters as a function of the other two, which can be resolved to 
eliminate one of them and reduce the new set to a five-parameter 
independent set equivalent to the original set. Hence, any trio 
of positive real numbers (a,b,c) satisfying (5.25) correspond to 
a 3rd-order loop design having the loop gain (K) and loop noise 
bandwidth (B n ) yet to be specified. 

Another equivalent set of design parameters applicable to 
3rd-order loops shall now be discussed. Its motivation arises from 
practical considerations: some 3rd-order loop realizations pre- 

sently employed in existing tracking systems exhibit a design 
formulation which can be identified as a special case of this 
equivalent set. The loop design parameters are now transformed 
using 




(5.26) 


which may be compared to (5.24) for reference purposes. The trio 
(x,y,z) is a dependent set that must be reduced to an independent 
pair, which along with (K, w 0 , t 2 ) yields the five-parameter de- 
sign set. For every trio (a,b,c) there exists a trio (x,y,z), 
which establishes the correspondence between B n and oj o (and between 
the two sets) from the relation 


0) 


<k) . 


/ x + — \ - - ( 1 . *1 _ Ii^o\ 

\ K ) y V y yK ) 




(5.27a) 
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*y - * a - j£> 

xy - z 


for Kx >> ai, 

0 

Ky >> t 2 uj 1 2 0 


(5.27b) 


The correspondence between the sets (K, B n , t 2 , a, b) and 
(K, o) o , x , x, y) for 3rd-order loops is analogous to that of the 
sets (K, B n , a) and (K, oj r ,c) for 2nd-order loops. The sets con- 
taining the (w 0 ) or (u) r ) parameters are based on critical (break) 
frequencies appearing in the closed-loop transfer function, while 
the sets containing the (B n ) parameter are based on additive noise 
locking error content through the loop noise bandwidth. All 
loops (perhaps of different order) with the same (B n ) value will 
exhibit the same rras additive noise error in the locked state, 
so that their relative merits can be established from dynamic lag 
error performance considerations.. The remaining loop parameters 
besides (B n ) can be selected to minimize the dynamic error with the 
knowledge that their selection will not alter the additive noise 
error content (which conversely happens with the other set). 

The steady-state dynamic lag error is again ft/K rad for a 
constant doppler shift (ft rad/sec) excitation, and infinite for 
a constant doppler rate (ft rad/sec 2 ) or higher dynamic inputs. 

These results are a consequence of the non-ideal integrators 
involved, since both steady-state errors above are zero for the 
ideal integrator case. The error transfer function is given by 


1 - H (s) 


l s + li S 2 + S 3 
K K K 


1 + (t, + i )s + (t j + -p-)s 2 + I*- s 

K K. K 


(5.28) 
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and the (s/K) and (t 2 s 2 /K) numerator terms are absent in the ideal 
integrator case. The steady-state error pattern may again be 
analyzed by long-division to yield the following results for the 
two input excitations in question: 


<t> (t) - - rad 

** K 


(5.29) 


and 




K 

!( 


t + T, - T 


t + T 




(5.30) 


The doppler rate error again increases linearly with time , 
so that a limit on the doppler rate duration or the tracking time 
is required to maintain lock. For small tracking times, the error 
may be approximated by its constant term and the loop analysis 
can proceed. Notice that the constraint t << t 2 valid for 2nd- 
order loops (see (5.20)) is not a-priori evident now. Moreover, 
even if such condition is assumed to hold, the constant term 
(i2x 2 /K) cannot be claimed to have any specific functional dependence 
on the loop noise bandwidth (B n ), since the time constant (t 2 ) 
represents a non- ideal integrator parameter in the loop filter. 

This implies that more specific information concerning the non- ideal 
integrator approximation in the loop filter realization should be 
included in the formulation to proceed with the 3rd-order loop 
analytical design. 
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5.3 The Mallinckrodt Third-Order Loo 


The generalized 3rd-order loop formulation presented in the 
previous section can now be used to account for the scarce 3rd-order 
loop realizations that have received significant attention in the 
literature. These particular designs have been motivated by cer- 
tain specific features they exhibit, but the aforesaid formulation 
provides a unified analytical approach from which these loops 

t- 

emerge merely as special cases of a wide class of 3rd-order loop 
realizations. Indeed, these few existing 3rd-order loop designs 
represent specific assignments for the 3rd-order loop parameters 
(a,b,c) or (x,y,z) previously discussed. The fact that these 
loops have been identified as possible solutions of the quadratic 
equation (5.25) considerably extends the number of potential 
3rd-order loop realizations and puts the few existing cases in the 
proper perspective. 

The main effort in 3rd-order loop analysis and design has 
concentrated on the Wiener, Pool and Mallinckrodt loops. The 
Wiener loop was derived from the analytical minimization of the 
total mean-square locking error (dynamic transient plus additive 
noise) in the presence of the doppler rate excitation plus white 
noise and using the loop linear model. The optimum loop filter 
and closed-loop transfer functions thus obtained yield a 3rd-order 

loop with two ideal integrators in the loop filter. The resultant 

( 3 9 27 \ 

— , ~, 1 or 

(x,y,z) = (2,2,1). This loop has received considerable attention 
yet the mean-square transient plus noise error can be severely 
questioned as a performance index as previously discussed. The 
interest in transient as opposed to steady-state errors is dubious 
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after acquisition has occurred and the loop maintains lock. Also, 
a bounded mean-square (time-averaged) transient error does not 
assure peak error overshoots below the locking values that validate 
the linear model. Finally, such transient mean-square error is not 
bounded when the non- ideal integrator approximations are acknowledged 
in practice, so that another performance index and optimization 
criterion becomes necessary. 

i 

The Pool and Mallinckrodt loop designs were motivated by 
certain dynamic transient response features (to the best of our 
knowledge) . They also exhibit two ideal integrators in their 
loop filters, thus sharing their dynamic tracking degradation 
when non-ideal integrators are employed in practice. Their loop 
designs correspond to the assignments (a,b,c) = (1/5.27, 1/6.91, 

1/6. 42) or (x,y,z) * (2 J~2~, 2, 1) in the Pool case, and (a,b,c) - 
(1/8.91, 1/19.85, 1/13.10) or (x,y,z) * (6,9,4) in the Mallinckrodt 
case. A comparison of the three 3rd-order loop designs in ques- 
tion is summarized in Table 2. Their basic distinction is reflected 
in the pole-zero patterns: the Wiener loop contaii. , complex conjugates 
in both transfer functions, the Pool loop only in H(s) but not in F(s), 
and the Mallinckrodt loop in neither. 

In any case, there exists a significant degradation in the 
dynamic error performance when the practical non-ideal integrators 
are acknowledged. The relative merits of any loop will then be 
governed by its effectiveness in reducing or maintaining such 
degradation within tolerable error bounds. This issue shall now 
be evaluated in detail for the Mallinckrodt loop realization only. 

This loup has been selected based on the following arguments: 
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TABLE 2 


Third-Order Loop Design with Ideal Integrators 


I 

V-n 

i 



Wiener 

Mallinckrodt 


1 + 2s ^ 2s 

1 + 6s + 9s’ 


w « <“0 

W 0 03 0 

H (s) 

i + — + ?£’+ ai 

1 + 5£ + 2il + 


(|J j. (1) Z to ’ 

0 0 0 

“o “1 


, , 2s , 2s z 

X + — + 

1 + «£ + 2ii 

KF(s) 

U). U) 1 

0 0 

"S 

s ’/«’ 

4s 2 /oj 3 


Bn 

5 

297 

w o 

3 

200 

(x.y,7.) 

(2.2.1) 

(6.9.4) 

(a.b.c) 

h . 9 . 27 \ 

/ i i 

\10 50 125/ 

\8.91 19.85 13 

pole-zero 
of H(s) 

2“.(» + f - i j 1 ) 

9-.(s + f-j 


(S + U ,)(s j^) 

4 ( S+ f )( S+ “ ’) 

pole-zero 
of KF (s ) 

2 tO o / , tO 0 + : w o\ 

s 2 ' 2 - 2 j 

^ (s + li)' 

4s z V 3 / 


4s 
u >? 


Pool 


1 + 2 72s + 2s ? 

,2 




<*) 


2 72s , 2s z . s 3 


I + 


(D, 


. .2 . , J 


1 + 2 + Zs^ 

w o 


s 2 /*i 


4 72 + 3 
4 72 - 1 

(2 72. 2. 1) 

(_L_ f _i_ , 1 \ 

\5.27’ 6.91 ’ 6.42/ 

2 4 + ^)’ 


(s + 0. 47 53a> 0 ) (s + 0. 7623w o ± j 

/ + ^ y 

s z V 72 / 


1 . 2376a) 0 ) 


I 


(a) It contains real poles and zeroes in F(s), so that the networks 
of Figure 5 are fully applicable; (b) It has been successfully 

'I 

developed and employed in existing tracking systems so that it 
represents state-of-the-art technology; (c) Its analytical develop- 
ment can be easily extended to the all-active case of Figure 5 (c) 
to reflect potential improvements that motivate its future imple- 
mentation. 

t 

The practical implemcr* ation of the Mallinckrodt loop shall 
be traced for itrf exist '..g realization based on the passive plus 
active case of Figure 5 (b) . The ideal vs. actual loop filter 
transfer functions are given by 


F.(s) - ( 1 + 3s/ui Q ) 

4K p K, s 2 /u>; (5.31) 


1 + R 2 C!S 1 + 

I 4 C 2 s 

1 + (Rj + R 2 )CxS 1 + 

R 3 (1 + K fl ) + R 4 J C 2 s 


(5.32a) 

1 + R 2 Ci s 1 + R 4 C 2 s 


1 ■+■ RjCjS 1 ■+■ K a R 3 C 2 s 

(5.32b) 


The loop filter zeroes are specified in terms of the critical 

frequency (w ) by the relation 
0 

R 2 Ci = R 4 c 2 = — (5.33) 

u o 

while the high frequency gain establishes the condition 



R, K 0 R 3 4K„K v 


(5.34) 
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which specifies the loop filter poles as 


pp _ p p / R 1 \ _ 3 / R i \ 3X 

R i c i R 2 L ’il ) l ) 

\R 2 / w . \ R 2 / 


(5.35a) 


K o R 3 C 2 


r 4 c 2 


/ K q R, \ . 3_ / K„rA _ 4K / R_2_ \ 

\ R4 / u 0 \ R 4 / 3^0 \ R 1 / 


4K 

3Xw! 


(5.35b) 


The doppler rate steady-state error of (5.30) can next be 
evaluated by noting that the loop filter realization corresponds to 


T 1 “ R 2 Cx + R4C2 “ — 


(5.36a) 


uj. 


t 2 = RiC x + K 0 R 3 C 2 - — + — - (5.36b) 

<*> 0 3Xuj J 

which yields the following results for the dynamic lag error contri- 
bution 


<|> ss = — rad (doppler shift) (5.37a) 

K 

<fc s = £ f t + — — 6 - + 1 ] rad (doppler rate)(5.37b) 

K \ u) Q 3 Xu)q K f 


The time dependent term in (5.37b) may be neglected under small 
tracking time conditions, and the constant term is noted to depend 
on the relative values of the (K, u) 0 , X) parameters. The conven- 
tional loop designs are based on large (K/ 0 ) v r .iues such that the 
(4K/3Xw*) term predominates, though the (3 X/co 0 ) term cannot be 
neglected in general. We shall assume that the doppler rate error 
can indeed be approximated by 4> ss = (4ft/3Xu)p rad and use X = 200 
to reflect existing loop designs. Under these conditions, the 
doppler rate error simplifies to <j> s 5 a (ft/150u)*) rad, which can be 
expressed in terms of the loop noise bandwidth u^ng B n = (297/200)w 0 
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for the Mallinckrodt loop. A comparison of this result with the 
2nd-order loop case shows a reduction in the 3rd-order loop dyna- 
mic error by a factor of 68 for the same noise bandwidth, or a 
reduction of the 3rd-order loop noise bandwidth by a factor of 
V68 for the same dynamic error. 
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6 . 


Range and Range Rate Data Errors 

The previous section has illustrated the various phase errors 
introduced in the range (tone) and range rate (carrier) signals by 
the received tracking loops. These signals are then fed to their 
respective range and range rate extractors that provide the raw 
tracking data samples corresponding to successive independent measure- 
ments. A phase shift A<J> on the reproduced tone induces a zero 

crossing (ranging event) time error At = — and range error 

u m 

Ar = j (At). Hence, any tone loop locking errors are directly con- 
verted into a range error via a (c/2w m ) multiplier. For example, 
the steady-state locking error to a tone doppler shift (0 m rad/sec) 
and doppler rate (ft m rad/ sec) were respectively derived as 

4> = ^ rad (6.1a) 

” K 

and 

Q 

<t> = 1J!L - rad (small tracking times) (6.1b) 

KB * 2 


where k = 1 for 2nd-order loops with C = 1/2, and < * 68 for 3rd- 
order Mallinckrodt loops with X = 200 and K/B n large. The corres- 
ponding dynamic range errors are then given by 


and 


H r 


c # 
2w m ' K 


r 

K 


(6.2a) 


Hr 


C . 

2(*) m kB^ 



(small tracking times) (6.2b) 


• • • 

where r and r represent the slant velocity (range rate) and accel- 
eration (range rate derivative) of the spacecraft. In turn, the 
additive noise locking error of (5.8) yields an rms range error 
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that can be expressed in terms of the tone loop SNR (regardless of 
the loop order) as 





2 /2w m 


(TNR)‘* 


(6.3) 


The range rate error is next evaluated by analyzing the two- 
way carrier doppler measurement operation. The doppler measure- 
ment can be reliably performed via cycle counting techniques: the 
number (N) of a bias (f^) plus doppler (f^) cycles in a time inter- 
val (T) satisfies the relation N = (f^ + f^ ) T, so that measurement 
of one parameter from the pair (N, T) while the other remains a 
known constant resolves the doppler effect as fj = - - f^. The ’bias 
frequency assures a positive value for (N) , if set larger than the 
peak doppler magnitude. The presence of doppler rates results in 
an average doppler measurement based on this procedure. The 
successive average doppler measurements (f^,^) then yield average 
range rate data r^ = — (f^,^) in discrete time samples (k) at a 
data rate of (1/T) samples per second. The choice of parameter to 
be measured and parameter to be fixed from the (N, T) pair yields 
the two range rate extractor realizations shown in Figures 6 and 7. 

The Fixed-N realization of Figure 6 measures the counting 
times (T^) in units (M^) quantized to the period (l/f r ) of a high 
frequency reference (f r ) . The first counter has a capacity of (N) 
cycle counts, and the 1st and (N + l)-th zero crossing of its bias 
plus doppler input trigger the start and stop pulses for the second 
counter. The number (M^) of reference cycles that occur in-between 
these two pulses then measures their time separation in (l/f r ) 
second units. The corresponding average doppler and range rate 
measurement are then obtained via the arithmetic computation 
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High Frequei y 
Reference Signal 



(a) Measurement Unit 



(b) Arithmetic Unit 


Figure 7 : Fixed-T Range Rate Extractor Realization 






(6.4) 





) 


The distinction between the measurement and arithmetic opera- 
tions is quite relevant for error analysis. The actual parameter 
being measured is M^, so that measurement errors can be accounted 
through their effects in altering this parameter. The presence of 
undesired phase jitter accompanying the bias plus doppler input 
shall advance or retard its zero crossings. The range rate measure- 
ment error can then be formulated from the reference error count 
AM, caused by a time shift in the 1st and (N + l)-th zero crossings, 
assuming the jitter is bounded so that cycle slippage effects can 
be neglected. Any trigger pulse jitter or reference cycle quanti- 
zation effects can be similarly analyzed to evaluate the measurement 
error induced. It should be emphasized that an equivalent phase 
jitter p(t) on the input signal alters the count only if there 
is a differential phase shift (A 4 >) t = <t> (T) - <t>(0) between the 1st 
and (N + l)-th zero crossing, since otherwise the latter are ad- 
vanced or retarded by the same amount. 

The reference count and range rate measurement errors intro- 
duced by a differential phase shift (Ap) T are given by 


and 


AM l = f (AT) = (A4>) t 

k r k 2*(f b + f d ) T 


“ ' 'if: < af d> * w ' *-<“> 


2f. 


h> t f d ( ) , _£ 

2 it N ' A,f 


4tt f c T 


(a*) t 


(6.5) 


( 6 . 6 ) 
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so that if the phase jitter 4>(t) i« assumed to be a WSS stochastic 
process with autocorrelation function R^(t), then the rms range 
rate error in a single measurement is given by 

9 ' -yfc [ R * <0) - R * (T) ] % (6 7a) 

a if %(T) « R*(0) (6.7b) 

The compact result of (6.7b) is obtained when the effective 
correlation time of the phase jitter process is small relative to the 
counting time (T) , so that the phase errors in the 1st and (N + l)-th 
zero crossings are independent and add mean-square wise. Such 
correlation time is governed by the shape of the autocorrelation 
function R^(t), which is affected by the equivalent phase filter 
processing the phase jitter. These filters correspond to the 
closed-loop transfer function of the previous section, so that the 
assumption of large B n T products validates (6.7b) and permits the 
following expression for additive white noise effects in terms of 
the carrier loop SNR when lock is maintained: 

a* » ^Tt < CNR V % (B " T lar 8 e) (6 8) 

The Fixed-T realization of Figure 7 measures the number (N' fc ) 

of bias plus doppler cycles that occur in a T-second interval. The 

first counter has start pulses every T-seconds, and registers the 

integer number (N^) of bias plus doppler cycles in this interval, 

but loses the fractional counts corresponding to the beginning (<5 ) 

k 

and end (p, ) of the interval as illustrated in Figure 8. The 
x 


- 63 - 


J 


i 


beginning time can be measured in real time using a second reference 

counter that starts and the beginning of the T-second interval and 

stops at the 1st zero crossing of the input signal. The number (M k ) 

of reference counts then measures this beginning time (6^) in (l/f r ) 

second units as 6^ = M k /f r • Conversely, the end time (y^) cannot be 

measured in real time since the first counter cannot tell a-priori 

which shall be the last (N k + 1) zero crossing to be registered. 

However, the assumption of a constant doppler effect permits to 

resolve the end time from the next beginning time as y * ^ - 6 , 

k f b +f d.k 

The net number of bias plus doppler cycles is thus given by 
NJ = N k + (f b + f d k ) (6^+ y^) , which can be expressed in terms of 
the measured parameters (N k , M fc , M k+] ) and equated to N k = (f fa + f rf k )T 
to resolve the doppler and range rate measurement as 






+ 1 


(M 


k+ > 



(6.9) 


It can again be noted that the actual measured parameters are counts 
(N k , M k , M k+) ) while the conversion to range rate data involves an 
arithmetic computation. The measurement errors can thus be accounted 
through their effects on these counts. The analysis is more compli- 
cated than in the previous case since now the zero crossings advanced 
or retarded by an input phase jitter can alter these counts in various 
cumulative or compensating ways. The possible cases and their error 
formulation are summarized in Figure 9 and Table 3, which conclude 
with the same range rate error expression (6.6) and corresponding 
rms value (6.7) as before. 

A basic distinction between the two cases lies in the parameter 
T appearing in these error expressions. The counting times (T k ) in 
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(k ♦ 1) T sec 


(k*?)T *«c 







Figure 9: Details of Fixed-T Measurement 
(AN* = 0 top and AN* f bottom) 



TABLE 3 



Fixed-T Measurement Error Cases 


Trans ition at ANj^ aM^ AM k+l 
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the Flxed-N case varies from measurement to measurement with the 
average doppler, so that's constant output data rate (1/T 0 ) can 
only be achieved by having T k < T„ and introducing dead zones of 
variable duration U k - T 0 - T k where the cycle counting is not 
performed. This effect is absent in the Fixed-T case, where T k - T # 
is applicable. The loss of available doppler information has 

motivated the terminology ’’destructive counting" for the Fixed-N 

, ** 

extractor, vs. "nondestructive counting" for the Fixed-T extractor. 
This dead zones usually eliminates error correlation over successive 
measurements in the Fixed-N case, while such correlation is present 
in the Fixed-T case and can be exploited in certain data smoothing 
procedures on the raw data samples obtained. 
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7 . Error Analysis and Performance Evaluatic * 

The formulation of the previous sections shall now be applied 
to derive quantitative range and range rate error results. However, 
the error expressions in question were contingent on the carrier 
and tone tracking loops maintaining lock, so that the loop track- 
ing performance must be evaluated along with the error magnitudes 
to validate the results. Also, the joint re uction of dynamic 
and noise locking error contributions involve a loop design compro- 
mise based on the loop noise bandwidth parameter. With reference 
to the carrier loop, the dynamic lag error of (6.1b) should be 
bounded to maintain lock in the presence of doppler rate tracking 
by allowing a sufficiently large loop bandwidth. The minimum 
bandwidth value that maintains lock should then be used to limit 
the rms noise locking error of (5.8) and the range rate error of 
(6.7b). The carrier loop design and range rate error formulation 
based on this approach is summarized in Table 4 for PM ranging. 

With reference to the tone loop design, the situation is more 
complicated because the dynamic lag locking error also contributes 
to range error. Two alternative design approaches shall be con- 
sidered : 

(a) The minimum loop bandwidth ir obtained by bound- 
ing the peak dynamic locking error of (6.1b), as in 
the carrier loop case. The loop bandwidth and dynamic 
locking error then become dependent on the tone fre- 
quency assignment through the doppler rate parameter 
(n„) . The dynamic range error of (6.2b) and the rms 
noise locking error of (5.8) will then also vary with 
the tone frequency since the loop bandwidth c'.oes. 
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The rms noise range error of (6.3) will exhibit a 
(oO" 3/4 functional dependence on the tone frequency 
rather than a simple dependence as a conse- 

quence of the loop bandwidth variation. The tone 
loop design and range error formulation based on 
this approach is summarized in Table 5 for PM ranging. 

(b) The minimum loop bandwidth is obtained by 
bounding the peak dynamic range error of (6.2b). The 
loop bandwidch and dynamic range error then become 
independent of the tone frequency assignment. How- 
ever, the peak dynamic locking error of (6.1b) will 
then vary with the tone frequency through the doppler 
rate parameter (ft m ) , unlike the previous case where 
the loop bandwidth dependence maintained such error 
invariant. The rms noise locking of (5.8) will also 
be independent of the tone frequency, while the rms 
noise range error of (6.3) will exhibit the simple 
(w m )"^ functional dependence. The tone loop design 
and range error f emulation based on this approach is 
summarized in Table 6 for PM ranging. 

The basic distinction between the two design approaches can 
be appreciated from the tabic notes. The dynamic locking error 
(y m ) is tone frequency independent in Table 5, but not in Table 6. 
Conversely, the noise locking error (o m ) is tone frequency depen- 
dent in Table 5, but not in Table 6. The design approach of Table 
5 may be supported from tracking considerations, since the loop 
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bandwidth is varied to account for different doppler rates to be 

tracked as different tone frequencies are used. The induced func- 

1/4 

tional dependence (w m ) of the noise locking error (o m ) then 
corresponds to the loop design compromise of wider bandwidths (more 
noise locking errors) if the same dynamic locking status is desired 
regardless of the tone frequency assignment. The design approach 
of Table 6 may be supported from range accuracy considerations, 
since it features the same dynamic range error (p r ) regardless of 
the tone frequency assignment, though at the expense of a variable 
dynamic locking status with the tone frequency. In both cases, the 
noise range error (a r ) decreases with higher tone frequencies but 
through a different functional dependence. The carrier locking 
error effects would contribute a (cos 2 A<J>) attenuation in the tone 
power level, and a (secA<j>) increase in the tone noise locking (o m ) 
and ranging (a,) errors. This effect is not included in the 
formulation for simplicity, so that the numerical results will be 
conditioned on maintaining good carrier lock, 

The tabulated formulation shall now be applied to perform a 
comprehensive range and range rate error analysis of the C0 2 laser 
ranging system. The link parameters under consideration have been 
specified in Table 7 for reference purposes. The carrier and tone 
sideband levels will be lower than the signal to noise density 
specified in the table by an amount governed by the modulation 
index (6) as previously discussed. The peak orbital acceleration 
shall be assumed to be r = 60.15 me "/sec 2 based on existing satellite 
(GEOS) orbital predictions. The peak carrier doppler rate is then 
h t = 224.56 tt x 10 5 rad/sec 2 for the laser frequency specified. The 
carrier loop and range rate error performance then follows from 
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the formulation of Table 4 upon assumption of peak dynamic locking 
error bounds. The three bounds given by y c = 1,5,10 degrees shall 
be assumed for both the 2nd (< = 1) and 3rd (< * 68) carrier loop 
realizations under consideration. An averaging time T = 1 second 
is assumed in the range rate extraction and rms additive noise error 
tabulation. 

The 2nd-order loop results are first presented in Table 8 
as a function of the modulation index. The tabulation of results 
was stopped when the loop CNR fell below 7 db, since the loop 
threshold is usually assumed in the 7-10 db range. It should be 
noted that the wider bandwidth (smaller dynamic locking error) case 
of Table 8(a) results in essentially below threshold operation even 
with all the signal power concentrated on the carrier component 
(no range data) . The situation is not significantly improved with 
the lower bandwidth (larger dynamic locking error) values of 
Tables 8(b) and 8(c), where only marginal above threshold levels and 
noise locking errors are obtained, and mainly with the signal power 
concentrated on the carrier component. The use of a narrower band- 
width to provide a reliable above-threshold performance is not 
acceptable, since the peak dynamic locking error would exceed re- 
liable locking bounds. These results thus force the conclusion 
that 2nd-order loops should not be employed for the carrier loop 
design. The 3rd-order loop results are next presented in Table 9 
following the same format. The CNR levels and corresponding noise 
locking errors are substantially improved to permit reliable above- 
threshold operation, except at the larger indeces where the signal 
power is concentrated on the tone at the carrier expense. 
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The tone loop performance and range error results are next 
presented. The dynamic locking error bound approach of Table 5 
yields the 2nd and 3rd order loop results of Tables 10 and 11, 
respectively. The bounds of y m = 1,5,10 degrees were again used, 
and the results were tabulated only for F m = 1 MHz for simplicity. 

In turn, the dynamic range error bound approach of Table 6 yields 
the 2nd and 3rd order loop results of Tables 12 and 13, respectively. 
The bounds of u r = 0.01 and 0.05 meters were used, and again the 
results were tabulated only for F m = 1 MHz for simplicity, The 
effect of varying the tone frequency is acknowledged in the notes 
as before, and the actual numerical error values obtained as the 
tone frequency varies shall be included in the performance curves 
that follow this tabulation. In all the cases it can be verified 
that the TNR is sufficiently high and the noise locking errors are 
sufficiently low to provide a reliable tone loop lock. 

The range rate error performance obtained with the 3rd-order 
carrier loop realization (Table 9) is illustrated in Figure 10 as 
a function of the modulation index (3) and the dynamic locking 
error bound (y c = 1,5,10 deg). The use of a weaker bound than 
10 deg is not recommendable , since there must be an error margin 
to account for additive and oscillator noise locking errors within 
the locked state. The use of a tighter bound than 1 deg is also 
undesirable, since the latter provides a sufficient error margin 
and a tighter bound would only increase the loop noise bandwidth, 
the additive noise locking error and the range rate additive noise 
error. The results illustrate that a low range rate error contri- 
bution from adcitive noise effects can be expected while maintaining 
a reliable dynamic error lock. The interest to account for 
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additional loop bandwidth constraints and potential carrier lock- 
ing and range rate error degradations from these results is further 
emphasized by the low error magnitude obtained. 

The range error performance obtained is next illustrated for 
both 2nd and 3rd order tone loop realizations, and also using both 
of the design approaches in Tables 5 and 6. The curves shall plot 
the rms additive noise error (a r ) vs. the modulation index (6) as 
a function of the tone frequency (F m ) employed. The relation 
y m (deg) = 2.4 F m (MHz) . y , (met) is always satisfied as previously 
noted, so that there is a one-to-one correspondence between the 
dynamic range error (y r ) and the tone frequency (F m ) in Table 5, and 
between the dynamic locking error (y m ) and the tone frequency (F m ) 
in Table 6. The tone frequency dependence of the plots then also 
becomes either a dynamic range error or a dynamic locking error 
dependence for the formulations in Tables 5 or 6 respectively. 

Hence, each curve in the (a r vs. 6) plots presented as a family 
corresponds to a pair (F m , y r ) or (F m , y m ) case. In turn, each 
family of curves presented in a different page corresponds to the 
option of 2nd vs. 3rd order loops (k = 1 or 68) and the option of 
a dynamic locking vs. ranging bound design approach (Table 5 or 6). 

The curves of Figures 11(a) - (f) are based on the formulation 
of Table 5, with the corresponding 2nd and 3rd order loop cases 
presented back to back to facilitate their comparison. The 
curves of Figures 12(a) - (d) are based on the formulation of 
Table 6, and follow the same format. It can be concluded that the 
peak dynamic and rms noise range errors (y r , a , ) can both be main- 
tained within the accuracy requirements using various moderate 
tone frequency and modulation index assignments, with either 2nd 
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or 3rd order tone loop realizations that maintain the peak dynamic 
and rms noise errors (|i„o n ) within locking bounds. The choice of 
modulation index should be also guided by the range rate error 
results of Figure 10. The range error curves are asymptotically 
flat at high indeces but rise steeply at low indeces, with the 
opposite pattern existing for the range rate error curves. Hence, 
a variation in the modulation index can be directly evaluated as 
an rms range vs. range rate error compromise, with the quantitative 
trade-offs given by the curves. Such trade-offs are a direct con- 
sequence of the power sharing compromise between the sidebands 
(tone) and carrier components of the PM ranging signal. 
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TABLE h 


Carrier Loop Design and Range Rate Error Formulation 


(1) Maximum acceleration: r = given (met/sec 2 ) 

• 2r 

(2) Maximum doppler rate: ft t = — oj, (rad/sec ) 


(3) 

Dynamic locking error bound: 

y — ^ w 

= chosen (rad) 

(A) 

Loop noise bandwidth: B„ =[ 

<B 2 

n 

ftc \ % = 

/2rw«\ 

k 

\ (Hz) 

(5) 

V 

Noise locking error: a = \ 

* V 

/N 0 B„ \ % 

\C<0) t j 

1 

r*— 

f 

/N.\ % / 2rw , ' 


‘ n/2 

\ / 

\2 

\P t / \C<U c J 

(6) 

Noise range rate error: = 

c 

*^2u) t T 

H 

O 

3 

II 

/ N 0 \ % / 2rw t ' 

\ P t / V C<u O 


Notes : 


(a) B„ directly proportional to 

(b) o t directly proportional to <D t 

V 

(c) o; directly proportional to w c for P„T large 

(d) CNR = for PM ranging 

N 0 B„ N 0 B„ 


k 

(rad rms) 


k 

(met/ sec) 
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TABLE 5 


Tone Loop Design and Range Error Formulation: Dynamic Locking Error Bound 


(1) Maximum acceleration: r = constant (met/sec 2 ) 

* 2r 

(2) Maximum doppler rate: (rad/sec 2 ) 


(3) Dynamic locking error bound 

(4) Loop noise bandwidth: B„ = 


y„ = = chosen (rad) 

<B 2 


/ D-. \ ^ _ / 2ru)„ \ 

\<P./ \c<y „ / 


(Hz) 


(5) Dynamic range error: y r = — — = (met) 

kB? 2w» 


(6) Noise locking error: a m 


(7) Noise range error: a, = 


2u)„ 


/2N.BA % 


i( 

N*. \* 

i 

3 

CM 

\ 4P, / 


2 \ 

P./ 

\CKy./ 

- ‘ — ( 

No. 

P, 

)* 

/ 2ru) „ \ 
\ck y m / 

% 

(met 


(rad rms) 


Notes: (a) y, = y„/( 2.4 F m ) (met) if U » (deg) and F. (MHz) 

(b) y „ independent of F„ 

(c) y , inversely proportional to F. 

(d) B„ directly proportional to F„^ 

(e) TNR inversely proportional to F„ ^ 

(f) o m directly proportional to F.^ 

(g) o . inversely proportional to F m 4 

(h) TNR = = 2J ‘ C B > P » 

NqBb N O B n 


for PM ranging (3 db processing gain) 


TABLE 6 


Tone Loop Design and Range Error Formulation: Dynamic Range Error Bound 


(1) Maximum acceleration: r = given (met/sec 2 ) 

* 2r 

(2) Maximum doppler rate: ft m = — a) m (rad/sec 2 ) 


(3) Dynamic range error bound: y, = = chosen (met) 

<B 2 

(4) Loop noise bandwidth: B n =(—^—\ (Hz) 


— 
<Ur / 


(5) Dynamic locking error: y m = 


. ft m _ 2y , a), 


(6) Noise locking error 


< B 2 c 

% 


(rad) 




(7) Noise range error: a, = 


2 to. 


a. 


\p./ V / 


(met rms) 


'Totes : (a) y m = 2.4 u f F B (deg) if y, (met) and F m (MHz') 

(b) y ,, B„, TNR, o m independent of F« 

(c) u „ directly proportional to F* 

(d) o r inversely proportional to F„ 


(e) TNR 


4P m 
2N e B„ 


1 (8)P. f or pm ranging (3 db processing gain) 


2N„ B 


I 


I 


TABLE 7 

Link Parameter Specifications 


Carrier Frequency 

2.8 x 10 1S GHz 

Transmitter Power 

23 dbw (200 w) 

Transmitter Antenna Gain 

100 db 

Receiver Antenna Gain 

101 db 

Transponder Reflector Gain 

162 db 

Path Loss (two-way) 

-493.78 db 

Miscellaneous Losses 

-26.40 db 

Receiver Signal Level 

-134.18 dbw 

Noise density 

-190 dbw - Hz 

Signal to Noise Density 

55.82 db/Hz 
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TABLE 8 

Range Rate Error Analysts (Second-Order Carrier Loop) 
8(a) Carrier Loop: FLL 2 , w c ■ 1 deg, B c - 63.92 KHz 


8 (rad 


CNR (db 





8(b) Carrier Loo 


CNR (db' 


.1 11.23 
.2 11.17 
.3 11.06 
.4 10.91 
.5 10.70 
.6 10.46 

1.7 10.16 

1.8 9.81 
'. 9 9.40 
.0 8.93 
.1 8.40 
.2 7.79 
.3 7.11 


6.61 

3.459 

x 10- 

6.74 

3.485 

x 10- 

6.95 

3.529 

x 10- 

7.26 

3.593 

x 10- 

7.66 

3.677 

i 

o 

H 

X 

5 deg, B t 

- 28.59 KHz 



o'. (met/sec 


11.11 

11.19 

11.34 

11.54 

11.81 

12.15 

12.58 

13.10 

13.72 

14.48 

15.40 

16.51 

17.87 
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I 


8(c) 

Carrier Loop: PLL 2 . u. 

- 10 deg, B t 

- 20.21 KHz 

( ra< *) 

CNR (db) 

Oc (deg) 

or (met/sec) 

0.1 

12.74 

9.343 

1.945 x 10“ 7 

0.2 

12.67 

9.413 

1.960 x 10" 7 

0.3 

12.56 

9.533 

1.985 x 10- 7 

0.4 

12.41 

9.704 

2.020 x 10- 7 

0.5 

12.21 

9.931 

2.068 x 10- 7 

0.6 

11.96 

10.22 

2.128 x 10“ 7 

0.7 

11.66 

10.58 

2.202 x 10- 7 

0.8 

11.31 

11.01 

2.293 x 10- 7 

0.9 

10.90 

11.54 

2.403 x 10- 7 

1.0 

10.44 

12.18 

2.536 x 10- 7 

1.1 

9.90 

12.95 

2.696 x 10- 7 

1.2 

9.30 

13.89 

2.891 x 10- 7 

1.3 

8.61 

15.03 

3.129 x 10- 7 

1.4 

7.83 

16.44 

3.423 x 10- 7 
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TABLE 9 

Range Rate Error Analysis (Third-Order Carrier Loop) 


9(a) 

Carrier Loop: PLLj, 

- 1 deg, B c 

- 7.751 KHz 

(rad) 

CNR (db) 

o c (deg) 

a; (met/sec) 

0.1 

16.90 

5.786 

1.205 x 10" 7 

0.2 

16.84 

5.829 

1.214 x 10 -7 

0.3 

16.73 

5.903 

1.229 x 10" 7 

0.4 

16.57 

6.009 

1.251 x 10" 7 

0.5 

16.37 

6.150 

1.280 x 10- 7 

0.6 

16.12 

6.328 

1.318 x 10“ 7 

0.7 

15.83 

6.549 

1.364 x lO -7 

0.8 

15.47 

6.819 

1.420 x lO -7 

0.9 

15.07 

7.147 

1.488 x 10“ 7 

1.0 

14.60 

7.542 

1.570 x lO -7 

1.1 

14.07 

8.020 

1.670 x 10 -7 

1.2 

13.46 

8.599 

1.790 x lO -7 

1.3 

12.77 

9.307 

1.938 x 10’ 7 

1.4 

11.99 

10.18 

2.120 x 10“ 7 

1.5 

11.11 

11.28 

2.348 x 10 -7 

1.6 

10.09 

12.67 

2.639 x 10‘ 7 

1.7 

8.92 

14.50 

3.019 x lO -7 

1.6 

7.55 

16.97 

3.534 x lO -7 
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9(b) 

Carrier Loop: PLLi; 

P ( - 5 deg, B { 

* 3.467 KHz 

( rad > 

CNR (db) 

< d *g) 

a: (met/sec) 

0.1 

20.40 

3.869 

0 8056 x 10‘ 7 

0.2 

20.33 

3.898 

0.8116 x 10“ 7 

0.3 

20.22 

3.948 

0.8219 x 10 -7 

0.4 

20.07 

4.019 

0.8367 x lO -7 

0.5 

19.87 

4.112 

0.8562 x 10 -7 

0.6 

19.62 

4.232 

0.8811 x 10' 7 

0.7 

19.32 

4.380 

J. 9119 x lO -7 

0.8 

18.97 

4.560 

0.9495 x 10‘ 7 

0.9 

18.56 

4.779 

0.9951 x lO -7 

1.0 

18.09 

5.044 

1.050 x lO" 7 

1.1 

17.56 

5.363 

1.117 x 10‘ 7 

1.2 

16.95 

5.751 

1.197 x 10- 7 

1.3 

16.27 

6.224 

1.296 x 10" 7 

1.4 

15.49 

6.809 

1.418 x 10" 7 

1.5 

14.60 

7.541 

1.570 x 10" 7 

1.6 

13.59 

8.475 

1.764 x 10" 7 

1.7 

12.42 

9.697 

2.019 x lO" 7 

1.8 

11.05 

11.35 

2.363 x 10' 7 
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2.451 KHz 


9(c) 

Carrier Loop: PLL 3 ; 

y t = 10 deg, B c 

= 2.451 KHz 

(rad) 

CNR (db) 

Or (deg) 

a: (met/ sec) 

0.1 

21.90 

3.254 

0.6774 x 10-' 

0.2 

21.84 

3.278 

0.6825 x 10- 

0.3 

21.73 

3.320 

0.6921 x 10- 

0.4 

21.57 

3.379 

0.7036 x 10" 

0.5 

21.37 

3.458 

0.7200 x 10- 

0.6 

21.12 

3.559 

0.7409 x 10' 

0.7 

20.83 

3.683 

0.7668 x 10' 

0.8 

20.47 

3.835 

0.7984 x 10' 

0.9 

20.07 

4.019 

0.8368 x 10' 

1.0 

19.60 

4.241 

0.8830 x 10' 

1.1 

19.07 

4.510 

0.9390 x 10- 

1.2 

18.46 

4.836 

1.007 x 10- 

1.3 

17.77 

5.234 

1.090 x 10- 

1.4 

16.99 

5.725 

1.192 x 10' 

1.5 

16.11 

6.341 

1.320 x 10* 

1.6 

15.09 

7.126 

1.484 x 10' 

1.7 

13.92 

8.155 

1.698 x 10' 

1.8 

12.55 

9.546 

1.987 x 10- 


C N - 2_ ' 


J 


TABLE 10 

Range Error Analysis (Second-Order Tone Loop, Dynamic Locking Error Bound) 


10(a) Tone Loop ; PLL 2 , y„ = 1 deg, F« = 1 MHz, B M = 12.02 Hz 


(rad) 

TNR (db) 

a* (deg) 

a, (met) 

0.1 

22.00 

3.218 

1 340 

0.2 

27.99 

1.615 

0.6725 

0.3 

31.45 

1.083 

0.4511 

0.4 

33.88 

0.8198 

0.3413 

0.5 

35.72 

0.6633 

0.2762 

0.6 

37.18 

0.5605 

0.2334 

0.7 

38.37 

0.4884 

0.2034 

0.8 

39.37 

0.4357 

0.1814 

0.9 

40.20 

0.3958 

0.1648 

1.0 

40.90 

0.3652 

0.1521 

1.1 

41.49 

0.3412 

0.1421 

1.2 

41.98 

0.3225 

0.1343 

1.3 

42.38 

0.3078 

0.1282 

1.4 

42.71 

0.2965 

0.1235 

1.5 

42.96 

0.2880 

0.1199 

1.6 

43.14 

0.2820 

0.1174 

1.7 

43.26 

0.2781 

0.1158 

1.8 

43.32 

0.2763 

0.1151 


Note: n, * y„/(2.4 F„) - 0.416 met for y„ = 1 deg and F„ * 1 MHz. 

1 ) 


F m scaling decreases (y, , TNR 2 , o,^ ) and increases (B 2 , 
proportionately, (y m remains invariant). 


a 




10(b) Tone Loop : PLL 2 y „ = 5 deg, F m = 1 MHz, B„ = 5.375 Hz 


(rad) 

TNR (db) 

(deg) 

a, (met) 

0.1 

25.49 

2.152 

0.8960 

0.2 

31.48 

1.080 

0.4497 

0.3 

34.95 

0.7245 

0.3017 

0.4 

37.37 

0.5482 

0.2283 

0.5 

39.21 

0.4436 

0.1847 

0.6 

40.67 

0.3748 

0.1561 

0.7 

41.87 

0.3266 

0.1360 

0.8 

42.86 

0.2914 

0.1213 

0.9 

43.69 

0.2647 

0.1102 

1.0 

44.39 

0.2442 

0.1017 

1.1 

44.98 

0.2282 

0.09502 

1.2 

45.47 

0.2157 

0.08980 

1.3 

45.88 

0.2059 

0.08571 

1.4 

46.20 

0.1983 

0.08256 

1.5 

46.46 

0.1926 

0.08020 

1.6 

46 . 64 

0.1886 

0.07851 

1.7 

46.76 

0.1860 

0.07745 

1.8 

46.82 

0.1848 

0.07695 

Mr = y m / (2. 

, 4F m ) = 2.08 met for 

y m = 5 deg and F„ 

= 1 MHz. 

F m scaling 

decreases (y r , TNR 2 

4A 

, o, ) and increases 


(B 2 ,o*) proportionately, (y m remains invariant). 
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10(c) Tone Loop : PLL 2 , y* 


10 deg, F 


1 MHz, B 


3.801 Hz 


(rad) 

TNR (db) 

o m (deg) 

a, (met) 

0.1 

27.00 

1.810 

0.7535 

0.2 

32.99 

0.9082 

0.3781 

0.3 

36.45 

0.6093 

0.2537 

0.4 

38.88 

0.4610 

0.1919 

0.5 

40.72 

0.3730 

0.1553 

0.6 

42.18 

0.3152 

0.1312 

0.7 

43.37 

0.2747 

0.1144 

0.8 

44.37 

0,2450 

0.1020 

0.9 

45.20 

0.2226 

0.09269 

1.0 

45.90 

0.2054 

0.08550 

1.1 

46.49 

0.1919 

0.07990 

1.2 

46.98 

0.1814 

0.07551 

1.3 

47.38 

0.1731 

0.07208 

1.4 

47.71 

0.1667 

0.06943 

1.5 

47.96 

0.1620 

0.06744 

1.6 

48.14 

0.1586 

0.06602 

1.7 

48.26 

0.1564 

0.06512 

1.8 

48.36 

0.1554 

0.06470 


Note: y, = y m /(2.4 F m ) = 4.16 met for y m « 10 deg and F m = 1 MHz. 

4 y 

F m scaling decreases (y r , TNR 2 , a, ) and increases 
(B 2 , oS) proportionately, (y m remains invariant). 
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TABLE 11 


Range Error Analysis (Third-Order Tone Loop, Dynamic Locking Error Bound) 
11(a) Tone Loop : PLLj, y „ * 3 deg, F„ - 1 MHz, B. - 1.458 Hz 


6 (rad) 

TNR (db) 

g - J de g) 

a. (met) 

0.1 

31.16 

1.121 

0.4666 

0.2 

37.15 

0.5624 

0.2342 

0.3 

40.62 

0.3773 

0.1571 

0.4 

43.04 

0.2855 

0.1189 

0.5 

44.88 

0.2310 

0.09618 

0.6 

46.34 

0.1952 

0.08127 

0.7 

47.54 

0.1701 

0.07082 

0.8 

48.53 

0.1517 

0.06317 

0.9 

49.36 

0.1378 

0.05740 

1.0 

50.06 

0.1272 

0.05295 

1.1 

50.65 

0.1188 

0.04948 

1.2 

51.14 

0.1123 

0.04676 

1.3 

51,55 

0.1072 

0.04463 

1.4 

51.87 

0.1033 

0.04299 

1.5 

52,12 

0.1003 

0.04176 

1.6 

52.31 

0.09819 

0.04089 

1.7 

52.43 

0.09686 

0.04033 

1.8 

52.43 

0.09623 

0.04007 

Note: y, * 

y„/(2.4 F. ) “ 0.416 met for y„ - 1 deg 

and F. * 1 MHz 

F„ scaling decreases (y, , 

TNR 2 , 0,^ ) and 

inert ases 

(B’ t 

a!) proportionally, 

(y . remains invariant). 


- 87 - 


11(b) Tone Loop : PLL 3 , y 


5 deg, F, 


1 MHz, B 


0.6518 Hz 


< rad ) 

TNR (db) 

Om (deg) 

a (met) 

0.1 

34.66 

0.7494 

0.3120 

0.2 

40.64 

0.3761 

0.1566 

0.3 

44.11 

0.2523 

0.1051 

0.4 

46.53 

0.1909 

0.07949 

0.5 

48.37 

0.1545 

0.06432 

0.6 

49.84 

0.1305 

0.05435 

0.7 

51.03 

0.1137 

0.04736 

0.8 

52.02 

0.1015 

0.04255 

0.9 

52.86 

0.09218 

0.03838 

1.0 

53.56 

0.08504 

0.03451 

1.1 

54.15 

0.07947 

0.03309 

1.2 

54.64 

0.07510 

0.03127 

1.3 

55.04 

0.07169 

0.02985 

1.4 

55.37 

C. 06905 

0.02875 

1.5 

55.62 

0.06707 

0.02793 

1.6 

55.80 

0.06567 

0.02734 

1.7 

55,92 

0.06477 

0.02697 

1.8 

55.98 

0.06435 

0.02680 


Note: u,«P./(2,4 F.) ■ 2.08 met for y. - 5 deg and F„ - 1 MHz. 
F„ scaling decreases (y, , TNR 2 , ) and Increases 

(B 2 , o!) proportionately, (V„ remains invariant). 
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11(c) Tone Loop : PLLs , M 


10 deg, F, 


1 MHz, B 


0.4609 Hz 


(rad) 

TNR (db) 

a. (deg) 

Or (met) 

0.1 

36.16 

0.6302 

0.2624 

0.2 

42.15 

0.3163 

0.1317 

0.3 

45.62 

0.2122 

0.08834 

0.4 

48.04 

0,1605 

0.06684 

0.5 

49.88 

0.1299 

0.05408 

0.6 

51,34 

0.1098 

0.04570 

0.7 

52.54 

0.09565 

0.03983 

0.8 

53.53 

0.08532 

0.03552 

0.9 

54.36 

0.07752 

0.03228 

1.0 

55.06 

0.07151 

0.02978 

1.1 

55,65 

0.06683 

0.02782 

1.2 

56.14 

0.06315 

0.02630 

1.3 

56.55 

0.06028 

0.02510 

1.4 

56, 8"" 

0.05807 

0.02418 

1.5 

57.12 

0.05640 

0.02348 

1.6 

57.31 

0.05522 

0.02299 

1.7 

57.43 

0.05447 

0.02268 

1.8 

57.48 

0.05411 

0.02253 


Note: Mr - Mm/ (2.4 F* ) * 4.16 met for M. « 10 deg and F. - 1 MHz. 
F. scaling decreases (M, , TNR 2 , o, /* ) and increases 
(B 2 , a!) proportionately, (M. remains invariant). 
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TABLE 12 


Range Error Analysis (Second-Order Tone Loop, Dynamic Range Error Bourn 


12(a) Tone Loop ; PLL 2 , y f ■ 0.01 met, B. - 77.56 Hz, F. - 1 MHz 


< rad > 

TNR (db) 

o m (deg) 

a, (met) 

0.1 

13.90 

8.174 

3.404 

0.2 

19.89 

4.102 

1.708 

0.3 

23.36 

2.752 

1.146 

0.4 

25.78 

2.082 

0.8671 

0.5 

27.62 

1.685 

0.7016 

0.6 

29.08 

1.424 

0.5928 

0.7 

30.28 

1.241 

0.5166 

0.8 

31.27 

1.107 

0.4608 

0.9 

32.10 

1.006 

0.4187 

1.0 

32.80 

0.9276 

0.3862 

1.1 

33.39 

0.8668 

0.3609 

1.2 

33.88 

0.8192 

0.3411 

1.3 

34.29 

0.7820 

0.3256 

1.4 

34.61 

0.7532 

0.3136 

1.5 

34.8b 

0.7316 

0.3046 

1.6 

35.05 

0.7163 

0.2982 

1.7 

35.17 

0.7065 

0.2942 

1.8 

35.22 

0.7020 

0.2923 


Note: * 2.4 y, F,, * 0.024 deg for y, * 0.01 met and F. * 1 MHz. 

F„ scaling increases y„ and decreases o, proportionately 
(y, , B„, TNR, <J„ remains invariant). 


- 90 - 


12(b) Tone Loop : PLL 2 , y 


0.05 met, B 


34.68 Hz, F 


1 MHz 


(rad) 

TNR (db) 

o„ (deg) 

a, (met) 

0.1 

17.40 

5.466 

2.276 

0.2 

23.38 

2.743 

1.142 

0.3 

26.85 

1.840 

0.7663 

0.4 

29.27 

1.393 

0.5798 

0.5 

31.11 

1.127 

0.4692 

0.6 

32.58 

0.9521 

0.3965 

0.7 

33.77 

0.8297 

0.3455 

0.8 

34.76 

0.7401 

0.3082 

0.9 

35.60 

0.6724 

0.2800 

1.0 

36.30 

0.6203 

0.2583 

1.1 

36.89 

0.5797 

0.2414 

1.2 

37.38 

0.5478 

0.2281 

1.3 

37.78 

0.5229 

0.2177 

1.4 

38.11 

0.5037 

0.2097 

1.5 

38.36 

0.4893 

0.2037 

1.6 

38.54 

0.4790 

0.1994 

1.7 

38.66 

0.4725 

0.1967 

1.8 

38.72 

0.4694 

0.1955 


Note: y„= 2.4 y r F„ = 0.12 deg for y r - 0.05 met and F„ * 1 MHz. 
F m scaling increases y,, and decreases a r proportionately 
(y r , B m , TNR, a „ remains invariant). 


- 91 - 


TABLE 13 


Range Error Analysis (Third-Order Tone Loop, Dynamic Range Error Bound) 


13(a) Tone Loop : PLLa, y r = 0.01 met, B* - 9.405 Hz, F„ - 1 MHz 


B 


Note : 


(rad) 

TNR (db) 

o m (dee) 

a r (met) 

0.1 

23.06 

2.847 

1.185 

0.2 

29.05 

1.429 

0.5948 

0.3 

35.52 

0.9584 

0.3991 

0.4 

34.94 

0.7252 

0.3019 

0.5 

36.78 

0.5867 

0.2443 

0.6 

38.24 

0.4958 

0.2064 

0.7 

39.44 

0.4321 

0.1799 

0.8 

40.43 

0.3854 

0.1605 

0.9 

41.26 

0.3502 

0.1458 

1.0 

41.96 

0.3230 

0.1345 

1.1 

42.55 

0.3019 

0.1257 

1.2 

43.04 

0.2853 

0.1188 

1.3 

43.45 

0.2723 

0.1134 

1.4 

43.77 

0.2623 

0.1092 

1.5 

44.03 

0.2548 

0.1061 

1.6 

44.21 

0.2494 

0.1039 

1.7 

44.33 

0.2466 

0.1024 

1.8 

44.39 

0.2444 

0.1018 

y - - 2.4 U r F m 

* 0.024 deg for 

y , * 0.01 met 

and F* - 1 HHz 


F» scaling increases y„ and decreases o r proportionately, 
(y r , B w , TNR, a m remains invariant). 


I 


J 


13(b) Tone Loop : PLL 3 , y, 85 0.05 met, B m * 4.206 Hz, F„ - 1 MHz 


(rad) 

TNR (db) 

o m (deg) 

a, (met) 

0.1 

26.56 

1.904 

0.7926 

0.2 

32. d5 

0.9554 

0.3978 

0.3 

36.01 

0.6409 

0.2669 

0.4 

38.44 

0.4849 

0.2019 

0.5 

40.27 

0.3924 

0.1634 

0.6 

41.74 

0.3316 

0.1381 

0.7 

42.93 

0.2889 

0.1203 

0.8 

43.93 

0.2577 

0.1073 

0.9 

44.76 

0.2342 

0.09750 

1.0 

45.46 

0.2160 

0.08995 

1.1 

'*■6 . 05 

0.2019 

0.08405 

1.2 

46 . 54 

0.1908 

0.07943 

1.3 

46.94 

0.1821 

0.07582 

1.4 

CM 

<* 

0.1754 

0.07304 

1.5 

47.52 

0.1704 

0,07094 

1.6 

47.70 

0.1668 

0.06945 

1.7 

47.82 

0.1645 

0.06851 

1.8 

47.88 

0.1635 

0.06807 


Note: Mm - 2.4 y, F m = 0.12 deg for u, * 0.05 met and F. - 1 MHz. 
F m scaling increases M m and decreases o. proportionately, 
(y f , B„, TNR, o. remains invariant). 
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Figure IT (a) Range Error Results (PLL - 2nd 
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Figure 12(d) Range Error Results (PLL = 3rd 
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Conclusions and Recommendations 


This report has covered in detail the work performed on the 
C0 2 laser ranging system. A system design based on independent 
range and range rate data was originally motivated, and the signal 
structure and subsystem processing were evaluated to this effect. 

The use of a high frequency range tone was recommended for the 
fine ranging operation, with lower subharmonic tones and a pseudo- 
noise code accounting for the ambiguity resolution operation. The 
fine ranging and ambiguity resolution operations should be per- 
formed sequentially to minimize the power sharing compromise be- 
tween the sidebands (range signal) and the carrier (range rate 
signal) components. 

The "ine range tone modulation on the carrier can be based on 
AM/DSB, P or coupling modulation methods, along with efficient 
coherent j. oduct demodulation techniques at the receiver. The 
carrier vs. sideband power sharing compromise was evaluated for 
these cases. The PM system features the advantages of a DSB narrow- 
band spectrum constant -envelope processing, plus an effective 
power sharing compromise. The use of SSB systems is undesired due 
to the harmful carrier locking error propagation effects. The 
coherent product demodulation exhibits phase-locked extraction of 
the carrier and demodulated tone components. The phase-locked loops 
act as lowpass filters on the input signal plus noise phase fluc- 
tuations, thus inducing a loop design compromise of wideband dyna- 
mic tracking potential vs. narrowband noise rejection capabilities. 

The loop design and tracking performance was evaluated for 2nd 
and 3rd order loop realizations. The dynamic tracking limitations 
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of non-ideal vs. ideal integrators in the loop filters were care- 
fully identified, and the locking error analysis accounted for the 
different nature of dynamic vs. noise error contributions. A con- 
ventional 2nd order loop and a state-of-the-art 3rd order loop 
were finally assumed as candidate realizations for the carrier 
and tone loop reproduction. The range and range rate accuracy 
limitations contributed by their dynamic tracking and additive 
noise errors were formulated in terms of the locking status and 
the SNR conditions prevailing in these loops. The loop design 
approach and resultant ranging error performance was then formulated 
to comply with both 2nd and 3rd order loop realizations, and to 
include the alternative control of dynamic locking vs. dynamic 
ranging error bounds . 

The rms range and range rate error contributions of additive 
noise were evaluated as a function of the modulation index, the 
tone frequency, the loop realization, and the dynamic locking or 
ranging error bound. The existing loop locking errors and SNR 
conditions were tabulated for lock verification purposes, and com- 
prehensive range and range rate error performance curves were pro- 
vided to guide the system design and performance compromises. A 
2nd-order carrier loop design could only provide below threshold 
or marginal above threshold locking conditions, and cannot be 
recommended. Conversely, a 3rd order carrier loop design exhibited 
a reliable locking performance, plus yielded a negligible range 
rate error contribution from additive noise effects . Both 2nd and 
3rd order tone loops were also capable of maintaining lock while 
yielding acceptable additive noise range errors for various moder- 
ate modulation index and tone frequency values. The logical 
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assignment of modulation parameters was then governed by the range 
vs. range rate accuracy compromise set by the system designer. 

The encouraging error performance thus obtained suggests the 
natural extension of the analysis to investigate the effects of 
reference and voltage-controlled oscillator noise. Both the 
carrier and tone frequencies involved in the laser ranging system 
are sufficiently high to motivate that such effects be accounted 
for in the locking and ranging error analysis. The oscillator 
phase noise inputs to the tracking loops will introduce additional 
locking constraints that may affect the resultant range and range 
rate error magnitudes . The oscillator noise errors will appear 
superimposed on the additive noise errors as an additive random 
disturbance, plus the additive noise errors may themselves vary due 
to the loop redesign that may be required to maintain a locked 
status . 

The carrier or tone loop tracking performance may also be 
improved by considering natural extensions of the 3rd-order loop 
realization assumed. The use of different parameter values in the 
Mallinckrodt loop filter approximation based on a passive plus 
active network should be investigated, along with the additional 
realistic improvements that can be expected with an all-active 
network realization. Even though such modifications have not been 
implemented (to the best of our knowledge) , they represent simple 
logical extensions of the state-of-the-art technology and merit 
investigation based on their theoretical locking and ranging error 
improvement potential. 
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9. Glossary of Basic Symbols 


A 

B c 

B m 

B„ 

FCs) 

F m 

H(s) 

•J. (6) 

K 

K. 

K P 

Kv 

L 

M 

M„ 

N 

N„ 

P 

Pc 

Pn 

P. 

T 

T 0 


ignal amplitude (volts) 
two-sided carrier loop noise bandwidth (Hz) 
two-sided tone loop noise bandwidth (Hz) 
two-sided loop noise bandwidth (Hz) 
loop filter transfer function 
tone frequency (MHz) 
closed- loop transfer function 
Bessel function coefficient 
loop gain (sec* 1 ) 
loop filter gain (volt/volt) 
phase detector gain (volt/rad) 
vco gain (rad/sec-volt) 
pseudo-noise code length (symbols) 
reference cycle count (cycles) 

2 

doppler dynamic magnitude (rad, rad/sec, rad/sec , etc.) 

doppler *»lus bias cycle count (cycles) 

one-sided additive thermal noise density (watts/Hz) 

pseudo-noise code period (sec) 

carrier component power (watts) 

tone component power (watts) 

total signal power (watts) 

doppler measurement interval (sec) 

output range rate data period (sec) 
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c - speed of light (met/sec, cm/sec) 

f - frequency unit (Hz) 

f b - bias frequency (Hz) 

f £ - carrier frequency (Hz) 

f d ■ doppler effect (Hz) 

f m - tone frequency (Hz) 

f, * reference frequency (Hz) 

r = slant range (met) 

r = slant velocity or range rate (met/ sec) 

•• j 

r = slant acceleration (met/ sec ) 
s = laplace operator 
t = time unit (sec) 
t p = propagation delay (sec) 

Q * doppler shift (rad/sec) 
n - doppler rate (rad/sec ) 

= carrier doppler shift (rad/sec) 

. 2 
fl c * carrier doppler rate (rad/sec ) 

= tone doppler shift (rad/sec) 

- tone doppler rate (rad/sec 2 ) 

a - AM sinusoidal modulation factor 
0 = PM sinusoidal modulation index (rad) 

y * coupling modulation factor 
k - 1 (2nd-order loop) or 68 (3rd-order loop) 
X » 3rd-order loop parameter 
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